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SUMMARY 

The  overall  objectives  of  this  work  are  to  (1)  refine  and  update  an  existing  soot  formation 
model  and  (2)  incoiporate  this  soot  model  into  a  code  describing  a  laminar,  opposed-jet 
diffusion  flame.  In  die  third  year  of  this  contract,  a  detailed  chemical  kinetic  mechanism 
for  the  pyrolysis  of  toluene  and  the  formation  of  polyaromatic  hydrocarbons  has  been 
developed.  This  chemical  kinetic  model  is  consistent  with  shock  tube  data  with  mass 
spectral  identification  of  intermediate  species.  These  data  were  obtained  in  the  first  year 
of  this  program.  This  kinetic  analysis  provides  a  better  understanding  of  how  aromatic 
rings  decompose  and  how  polyaromatic  species  grow.  In  a  related  effort,  the  previously 
developed  soot  formation  code,  based  on  a  sectional  aerosol  model  for  predicting  soot 
inception,  growth,  and  oxidation  in  a  premixed  flame  has  been  fully  integrated  into  an 
opposed-jet,  diffusion  flame  code.  The  new  code  includes  effects  due  to  radiation  from 
both  gaseous  species  and  particulates  as  well  as  scavenging  of  species  by  soot.  The  code 
treats  particle  transport  including  thermophoresis.  Predictions  from  a  low  strain  rate, 
lightly-sooting,  methane-fueled,  opposed-jet,  diffusion  flame  are  included  in  this  report. 
Perturbation  studies  demonstrate  the  importance  of  fully  integrating  soot  production, 
radiation,  and  scavenging  in  order  to  reasonably  predict  bulk  parameters  such  as 
temperature  as  well  as  species  concentrations,  sooting  levels,  and  radiation  loads. 


I.  Introduction 

The  overall  goal  of  this  three  year  effort  is  to  incorporate  a  soot  model  into  a  code  for  a 
laminar,  opposed-jet,  diffusion  flame.  Specific  objectives  necessary  to  reach  this  goal 
are:  (1)  identify  reduced  kinetic  mechanisms  which  adequately  describe  ring  formation 
and  growth  processes,  (2)  obtain  and  interpret  new  data  on  polyaromatic  hydrocarbons 
(PAH)  formed  during  hydrocarbon  pyrolysis,  (3)  refine  and  update  an  existing  soot 
formation  model,  (4)  incorporate  a  radiation  model  into  a  code  describing  a  laminar, 
opposed-jet  diffusion  flame,  and  (5)  incorporate  the  soot  model  into  the  flame  code.  In 
the  annual  reports  for  die  first  two  years  of  this  program  [1,2],  results  and  progress 
related  to  much  of  the  above  were  discussed.  In  particular,  (1)  a  sinqilified  kinetic 
sequence  describing  the  dominant  paths  for  benzene  formation  and  destruction  was 
developed  and  verified  by  comparison  to  experiment,  (2)  new  experimental  data  on 
toluene  pyrolysis  were  obtained  with  a  focus  on  die  identification  of  key  intomediate 
species  via  mass  spectrometry,  (3)  temperature-dependent  thermodynamic  parameters 
were  determined  for  many  intermediate  hydrocarbon  species  and  those  for  PAHs  were 
obtained  from  the  literature  or  calculated  using  known  correlations,  (4)  the  soot 
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formation  code  for  a  premixed  flame  was  modifled  to  include  effects  due  to  soot  aging, 
(5)  a  band  radiation  model  was  developed  to  treat  radiation  from  flames  with  varying 
degrees  of  optical  thickness,  (6)  predictions  from  the  sectional  model  for  soot  formation 
were  compared  with  predictions  from  global  soot  formation  models  and  limitations  of 
these  global  models  were  identified,  and  (7)  die  sectional  soot  formation  model  was 
added  to  the  opposed  jet  diffusion  flame  code,  although  debugging  of  the  combined 
codes  was  incomplete. 

In  this  final  report,  these  efforts  are  summarized  by  a  manuscript  describing  a  detailed 
chemical  kinetic  mechanism  for  toluene  pyrolytic  decay  and  PAH  formation  and  a  paper 
specifying  the  combined  soot  formation  and  diffusion  flame  codes.  These  manuscripts 
are  mostly  self-explanatory  so  only  some  of  the  highlights  are  reviewed  here.  In  the  case 
of  the  toluene  paper,  however,  some  details  on  the  chemical  kinetic  mechanism  and 
thermodynamic  data  could  not  be  included  in  die  publication  due  to  space  limitations; 
consequently,  this  infonnation  is  included  here  for  completeness. 


II.  Results 


II.  A.  Kinetics  of  Toluene  Pyrolysis 

Reaction  mechanisms  and  kinetics  for  toluene  decomposition  and  PAH  formation  are 
proposed  based  on  literature  results  and  new  experimental  data  obtained  from  the  single¬ 
pulse  shock  tube  (SPST).  These  results  are  described  in  Appendix  A  of  this  report.  The 
detailed  reaction  mechanism  is  reproduced  in  Table  1  and  die  thermodynamic  data 
required  for  the  kinetic  estimations  and  calculations  of  reverse  reaction  rates  are  provided 
in  Table  2. 

Key  results  from  this  study  are  the  identification  of  a  low  energy  path  for  die 
decomposition  of  benzyl  radical  based  on  product  distribution,  dieoretical  arguments,  and 
evidence  obtained  by  shock  heating  similar  chemical  systems,  specifically 
cyclopentadiene,  norbomadiene,  and  mixtures  of  cyclopentadiene  and  acetylene.  The 
stability  of  benzyl,  an  aromatic  radical,  is  resonandy  enhanced  and  the  products  of  its 
decomposition  have  eluded  researchers  for  many  years.  Evidence  for  and  specifics  of  this 
new  proposal  are  reviewed  here. 

Two  low  temperature  products  observed  during  toluene  pyrolysis  coincident  with  early 
formation  of  acetylene  (presumably  a  decomposition  product  of  benzyl  or  phenyl)  are 
cyclopentadiene  (cpd)  and  benzyl-cyclopentadienyl  (b-cpd).  The  latter  product  is  almost 
entirely  formed  during  the  quenching  process  via  the  recombination  of  the  two  stable 

radicals,  benzyl  and  cyclopentadienyl.  Traces  of  several  odier  C,  species  at  masses  66 
and  68  were  observed  at  levels  typically  5  to  10  times  below  those  of  cpd.  The  presence 
of  die  Cj  conqxiunds  supports  arguments  for  decomposition  of  benzyl  radical  into  a  - 
iqiecies.  Cpd  is  preferred  since  cyclopentadienyl  offers  the  lowest  energy  path  for 
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TABLE  I 

PROPOSED  MECHANISM  FOR  TOLUENE  PYROLYSIS 
log  k  =  log  A-  E/R/T/2J03,  units:  cc-  mole-  sec-  kcal 


Reactions 

Forward  Rate  Constant 

Reverse  Rate  Constant 

Log  A 

E 

Log  A 

E 

1  C7Hg=ben2yl+H 

m 

88.9 

■fM 

-0.2 

2  C7H8=phenyl+CH3 

97.0 

-0.8 

3  C7H8+H=benzyI+H2 

14.10 

8.4 

23.7 

4  C7H8+CH3=CH4+benzyl 

11.1 

13.08 

27.1 

5  CgH6=phenyl+H 

116.0 

14.46 

6.4 

6  C5H6+H=phenyl+H2 

14.40 

16.0 

12.71 

10.7 

7  C6H6+CH3=phenyl+CH4 

13.30 

16.0 

13.10 

11.5 

8  C7H8+H=C6H6+CH3 

13.08 

5.1 

12.30 

17.0 

9  2benzyl=bibenzyl 

12.70 

0.5 

14.93 

60.0 

10  H+biben2yl->stilbene+H2+H 

13.70 

5.0 

0.00 

0.0 

1 1  H+biben2yl->C8H8+phenyl+H2 

13.70 

13.0 

0.00 

0.0 

12  stilbene+H— »phenanthrene+H2+H 

14.40 

16.0 

0.00 

0.0 

13  stilbeneK)H3->phenanthrene+CH4+H 

13.30 

16.0 

0.00 

0.0 

14  stilbeQe=phenanthtene+H2 

8.00 

40.0 

8.37 

45.5 

15  2phenyl=C  12^10 

12.50 

0.0 

15.96 

107.5 

16phenyl-K:6H6=Ci2Hio+H 

12.00 

7.0 

13.16 

4.8 

17  C2H2+phenyl=C6H5C2H+H 

12.00 

5.0 

13.81 

6.6 

18  Ci2Hio+H=Ci2H9+H2 

14.40 

16.0 

13.04 

10.5 

12.30 

5.0 

16.42 

55.8 

20  phenyl+C7H8=DPM+H 

11.70 

7.0 

11.80 

8.2 

21  DPM+H->fluo^ene^-H+H2 

14.40 

16.0 

0.00 

0.0 

22  DPM=fhiorene+H2 

8.00 

32.0 

10.24 

18.3 

23  fluorene+CH3-»anthracene+H2+H 

12.30 

17.0 

0.00 

0.0 

24  beiizyl+C7H8=DMDP+H 

12.48 

28.0 

12.19 

3.3 

25  benzyl-K:7H8=DMDPl+H 

12.30 

28.0 

12.02 

3.3 

26  DMDPl+H-»dhanthracene^-H+H2 

14.00 

7.0 

0.00 

0.0 

27  DMDPl=dlianthracene+-H2 

8.00 

40.0 

9.67 

10.7 

28  dbaiithracene=antiiracene+H2 

12.36 

43.8 

12.02 

52.2 

29  aiitiiracene°^henantlirene 

12.90 

65.0 

12.73 

67.0 

30  phenanthrene+H=phenanthrenyl+H2 

14.40 

16.0 

13.45 

11.3 

31  phenanthrenyl+C2H2“*pyrene+H 

12.00 

5.0 

0.00 

0.0 

32  plienyl+C5H5C2H-^pheiianthrene+H 

12.00 

5.0 

0.00 

0.0 
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TABLE  1  (coatiBued) 

PROPOSED  MECHANISM  FOR  TOLUENE  PYROLYSIS 
log  k  log  A-  E/R/T/2  J03,  uaits:  cc-  mole-  sec-  kcal 


ReactioBs 

Forward  Rate  CoBStaBt 

Reverse  Rate  Coastaat 

Log  A 

E 

Log  A 

E 

33  C^H4C2H+C^H^C2H-^pyreae+H 

11.00 

5.0 

■■ 

0.0 

34  H+C7Hg=CH3C6H4+H2 

14.00 

16.0 

9.7 

35  CH3C6H4+beiizyl=DMDP 

12.18 

0.0 

mm 

86.0 

36  H+DMDP=DPM+CH3 

13.08 

5.1 

11.60 

22.4 

37  CH3C6H4+ben2yl=DMDPl 

12.18 

0.0 

14.29 

86.0 

38  CH3C6H4+C7H8=DMDP+H 

11.30 

7.0 

12.20 

3.9 

39  C7H8+CH3CH3C6H4+CH4 

13.00 

16.0 

12.70 

10.4 

40  CH3C6H4+C7H8=DMDP1+H 

11.18 

7.0 

12.08 

3.9 

41  CH3C^+C7H8=benzyl+C7H8 

11.48 

3.0 

12.66 

24.6 

42  CH3C5H4+C2H2-^indene+'H 

12.00 

5.0 

0.00 

0.0 

43  indene+CH3->naphthalene+-H2+H 

12.30 

17.00 

0.00 

0.0 

44  phenyl+C7H8=C6H6+benzyl 

12.20 

11.1 

13.29 

31.7 

45  CH3+C7H8=*xylene+H 

11.31 

17.0 

12.91 

4.3 

46  CH3+CH3C5H4*xylene 

13.30 

0.0 

17.30 

98.0 

47  benzyl+CH3=C8HjQ 

12.70 

0.0 

15.00 

73.4 

48  b€nzyl=c-C5H5+C2H2 

13.78 

70.0 

11.57 

8.3 

49  benzyI=C4H4+C3H3 

14.30 

83.6 

8.28 

-15.6 

50  phenyl=C4H3+C2*^2 

13.65 

72.5 

7.79 

-25.8 

51  CH3C6H4=C4H3+C3H4 

13.65 

72.5 

8.38 

-22.1 

52  C3H3+b«nzyl=naphthalene+^2H 

11.78 

0.0 

0.0 

0.0 

53  beiizyl+c-C5H5=b-cpd 

12.70 

0.0 

13.81 

49.3 

54  benzyl+C5H^=b-q)d+H 

11.70 

30.0 

11.87 

-0.3 

55  C7H8+c-C5H5=b-cpd+H 

11.70 

40.0 

11.59 

0.2 

56H+C8Hio=C6H6+C2H5 

13.08 

5.1 

12.14 

20.1 

57  H+C8Hio->C8H8+H2+H 

14.40 

8.4 

0.00 

0.0 

58  H+C8H8->C6H5C2H+H2+H 

14.84 

14.5 

0.00 

0.0 

59  C6H5C2H+H=C6H4C2H+H2 

14.40 

16.0 

13.18 

10.7 

60  C6H4C2H+C2H2=naphthalenyl 

12.00 

5.0 

16.54 

96.2 

61  naplithalene+H=iK^)hthalenyl+H2 

14.40 

16.0 

12.77 

10.5 

62  n^hthalenyl+C2H2~>8ceiiaphthylene+H 

12.48 

5.0 

0.00 

0.0 

63  iu^hthalenyl+C2H2=CioH7C2H+H 

12.00 

5.0 

13.81 

6.6 
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TABLE  1  (coBtiBued) 

PROPOSED  MECHANISM  FOR  TOLUENE  PYROLYSIS 
log  k  >  log  A-  E/R/T/2 J03,  uBits:  cc-  mole-  sec-  Real 


Forward  Rate  CoBstaot  Reverse  Rate  CoBstaot 


ReactioBS 


64  CioH7C2H+H=CioH6C2H+H2 

65  CioH5C2H+C2H2=phenanthrenyl 

66  naphthalenyl+CgHg— »pyrene+H+H2 

67  iiaphthalene+phenyl-4pyrene+H+H2 

68  H-K:8H8=C6H6+C2H3 

69  2CH3=C2H6 

70  C2H3+C7H8=C2H4+C7H7 

71  C3H5-K:2H2=C5H7 

72  C5H7C5H6+H 

73  C5H6=H-K;-C5H5 

74  2c-C5H5=bicpdyl 

75  H+C5H6=1'C5H7 

76  1-C5H7+C5H6=C5H8+c.C5H5 

77  C5H8+HC2H4+C3H5 

78  H+C5H6=H2+c-C5H5 

79  H+C5H6=C2H2+C3H5 
8OH+C5H6C2H4+C3H3 

81  c-C5H5=C2H2+C3H3 

82  2C3H3— >phenyl+H 

83  CH3+C5H5CH4+C-C5H5 

84  C3H5+C7H8=C3Hg+benzyl 

85  CH3CHCH+C7H8=C3H6+benzyl 
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TABLE  2 

THERMODYNAMIC  PROPERTIES  OF  SELECTED  SPECIES 

Name 

AH®, 

S® 

Cp 

(cal/mole/K) 

kcal/mole 

cal/mole/K 

300K 

500K 

800K 

lOOOK 

1500K 

Toluene 

12.0 

76.7 

24.8 

40.6 

55.9 

62.4 

71.6 

Benzyl 

47.8 

73.5 

26.2 

40.8 

54.5 

60.2 

68.2 

Mediylphenyl 

71.4 

81.7 

24.6 

38.7 

52.6 

58.4 

67.0 

Ethylbenzene 

70.7 

86.1 

30.8 

49.3 

67.1 

74.8 

85.9 

Styrene 

35.2 

82.5 

29.4 

45.9 

60.9 

67.8 

77.5 

Xylene 

43.6 

84.3 

30.2 

48.5 

67.3 

75.5 

84.8 

Bibenzyl 

34.3 

113.5 

49.3 

79.8 

108.1 

120.1 

136.6 

Trans-stilbene 

56.0 

108.3 

48.0 

73.7 

102.0 

113.6 

128.1 

Dimethyldiphenyl* 

30.5 

121.3 

48.0 

77.9 

106.8 

119.4 

137.4 

Diphenylmethane** 

33.0 

112.3 

42.4 

70.3 

96.5 

107.6 

123.5 

Benzylcyclopentadienyl 

52.4 

101.8 

41.7 

70.2 

95.1 

105.1 

119.6 

Biphenyl 

43.6 

93.7 

39.7 

64.4 

87.0 

96.0 

108.2 

Naphthalene 

36.0 

79.5 

31.8 

52.0 

70.6 

78.6 

90.0 

Indene 

39.1 

80.7 

29.6 

48.2 

66.1 

73.4 

84.0 

Fluorene 

47.0 

91.4 

40.1 

67.6 

92.2 

101.8 

115.0 

Acenaphthylene 

61.7 

67.2 

38.0 

60.7 

81.0 

89.1 

100.7 

Dihydroanthracene 

58.4 

100.6 

47.3 

76.7 

103.9 

115.1 

130.6 

Anthracene 

52.1 

94.0 

44.4 

72.4 

97.3 

107.2 

121.3 

Phenanthrene 

50.0 

94.5 

44.9 

72.7 

97.3 

107.2 

121.3 

Pyrene 

55.2 

96.5 

48.9 

79.7 

106.7 

117.4 

132.3 

*  DMDP  and  DMDPl 
♦*DPM 
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decomposition  (64  kcal/mole).  The  mechanism  for  this  process  has  not  been  described 
wel].  Thus,  we  explored  the  reverse  process  by  comparing  results  from  cpd  pyrolysis 

with  those  of  cpd/C2H2  pyrolysis.  The  rate  of  toluene  formation  was  enhanced  by  about 
an  order  of  magnitude  in  the  presence  of  acetylene.  Following  a  suggestion  by  Kiefer,  the 
toluene  fomoation  was  then  compared  with  the  reverse  of  the  known  elimination  of 
acetylene  from  norboraadiene  (nbde).  Nbde  should  be  an  intermediate  in  the  conversion 
between  toluene  and  cpd  plus  acetylene: 

c-CjH^  +  C^Hj  <=>  nbde  <=>  toluene. 

Rate  constants  for  the  formation  of  toluene  by  the  above  process  were  calculated  from 
the  experimental  data  and  were  compared  to  literature  values  for  the  reverse  process. 
After  reducing  the  reverse  rate  by  a  factor  of  four  to  account  for  decomposition  of  nbde 
into  the  reactants  and  isomerization  to  cycloheptatriene,  the  extrapolated  rate  for  the 
addition  reaction  is  more  than  a  factor  of  two  below  the  experimental  val'ie.  More 
importantly,  the  temperature  dependence  of  the  rate  constant  of  the  published  value  for 
the  reverse  process  (25  kcal/mole)  is  significantly  lower  than  our  experimental  value  of 
nearly  60  kcal/mole.  Alternatively,  we  may  be  observing  the  radical 
addition/isomerization: 

c-CjHj  +  CjHj  <=»  norbomadienyl  <=>  benzyl. 

The  possibility  of  the  reverse  of  this  process  has  been  suggested  previously.  Thus,  the 
observed  temperature  dependence  would  be  due  principally  to  the  temperature 
dependence  of  cpd  decomposition  and  formation  of  cyclopentadienyl.  The  activation 
energy  for  the  deconqiosition  of  benzyl  into  cyclopentadienyl  as  listed  in  Table  1  (70 
kcal/mole)  is  based  on  low  barriers  for  the  brit^ging  process  to  form  norbomadienyl  and 
subsequent  acetylene  elimination.  During  the  toluene  experiments,  a  trace  peak  (at  mass 
92)  was  observed  which  eluted  near  the  retention  time  for  nbde,  but  insufficient  mass 
was  recovered  for  positive  identification.  To  provide  some  further  confirmation  of  these 
types  of  processes,  nbde  was  pyrolyzed  in  the  SPST.  Even  at  the  lowest  temperatures 
tested,  near  1  lOOK,  only  traces  of  nbde  remained  after  pyrolysis,  while  the  products  cpd, 
toluene,  cycloheptatriene,  and  acetylene  were  observed.  The  rapid  loss  of  nbde  was 
consistent  with  known  rates  for  nbde  isomerization/decomposition.  It  is  difficult  to 
extract  information  from  this  latter  study  with  nbde  to  confirm  the  proposed  benzyl 
radical  decomposition  pathway. 

Nevertheless,  Reaction  48  was  adopted  as  the  preferred  low  temperature  benzyl 
decomposition  path.  An  energy  diagram  for  this  process  is  provided  in  Fig.  5  of 
Appendix  A.  The  cpd  and  b-cpd  production  were  used  to  help  determine  a  rate  for  benzyl 
radical  decomposition.  Unfortunately,  a  sensitivity  analysis  indicates  that  the 
concentrations  of  these  species  are  dependent  on  both  the  rates  for  decomposition  of 
benzyl  (Reaction  48)  and  for  cyclopentadienyl  radicals  (Reaction  81).  The  rates 
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suggested  in  Table  1  for  these  two  reactions  should  be  considered  accurate  to  about  a 
factor  of  three  in  the  temperature  range  of  1 300  to  1 500K. 

Other  key  conclusions  from  the  toluene  paper  include  identification  of  several  ring 
growdi  processes.  A  key  stq)  not  well  recognized  previously  are  four-centered  reactions 

causing  ring  closure  and  direct  elimination.  Pathways  involving  sequential  acetylene 
addition  to  aromatic  species  are  not  fast  enough  to  explain  the  observed  formation  rates 
of  PAH.  In  addition,  aromatic  isomerization  reactions  were  demonstrated  with  half-lives 
at  1600K  near  70  microseconds. 


II.  B.  Soot  Formation  Modeling  in  Diffusion  Flames 

A  model  describing  soot  formation  in  an  opposed  jet  diffusion  flame  has  been  developed 
by  coupling  the  dynamical  sectional  equations  for  spheroid  growth  with  the  conservation 
equations  and  complex  chemistry  of  counterflow  diffusion  flames.  A  schematic  of  die 
model  indicating  the  fully  integrated  nature  of  the  soot  formation  code  and  the  flame 
model  is  provided  in  Fig.  1.  A  diagram  of  the  opposed  jet  flame  modeled  in  this  study  is 
shown  in  Fig.  2  and  details  of  the  model  are  described  in  Appendix  B  of  this  report.  The 
particle  kinetics  include  inception,  surface  growth  and  oxidation,  and  coalescence  within 
the  framework  of  the  discrete  size  spectrum,  or  sectional  representation.  Model 
predictions  demonstrate  that  the  particle  growth  dynamics  and  transport  (by 
thermophoresis,  diffusion,  and  convection),  gas  phase  chemistry,  flame  radiation,  and 
flow  hydrodynamics  are  all  strongly  coupled.  Particle  inception  as  well  as  surface  growth 
and  oxidation  are  all  calculated  from  local  conditions,  including  concentrations  of  key 
species.  As  part  of  this  modeling  effort,  the  benzene  formation/decomposition 
mechanism  was  updated  and  a  steady-state  solution  for  predicting  formation  rates  of 
polyaromatic  species  was  developed.  Scrubbing  of  gas  phase  species  during  growth  and 
oxidation  of  the  particles  is  included  in  the  analysis.  Non-adiabatic  radiative  loss  by  gas 
and  particle  emission  is  also  coupled  to  the  energy  eqtiation.  In  Figs.  3  and  4,  soot 
volume  fractions  are  plotted  against  the  normalized  coordinant,  z',  as  defined  in  Fig.  2 
and  (mathematically)  in  Appendix  B.  As  can  be  seen,  predictions  are  qualitatively  in 
good  agreement  with  experimental  data.  (For  details  of  the  calculations  and  die 
experimental  data  see  Appendix  B  and  Ref.  [3],  respectively).  Soot  appears  first  just  on 
the  fuel-rich  side  of  the  peak  temperature  region,  near  the  benzene/acetylene  peak,  and 
grows  to  a  maximum  near  the  stagnation  plane.  A  strong  coupling  between  thermal 
radiation  and  soot  formation  was  observed  since  gas  cooling  reduces  acetylene  and 
benzene  formation  as  well  as  the  Arrhenius  factors  in  inception  and  surface  growth. 
Under  high  surface  growth  conditions,  total  soot  production  becomes  less  sensitive  to 
soot  inception  rates  but  can  strongly  affect  local  concentrations  of  gaseous  species  such 
as  acetylene.  This  depletion  of  acetylene  will  in  turn  reduce  the  rate  of  surface  growth. 
Oxidation  can  significantly  affect  predictions  of  soot  growth.  The  algorithm  is  highly 
efficient  with  accurate  solutions  of  soot  volume  fractions  obtained  with  only  two  or  three 
size  classes,  although  accurate  solutions  of  average  size  and  number  density  require  15  to 
20  size  classes. 
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Fig.  1  Diagram  of  Integrated  Soot  Formation 
Model  and  Flame  Gode 
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g.  2  Configuration  of  Opposed  Jet,  Laminar  Diffusion  Flame 
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in.  List  of  Publications 

A  manuscript  entitled  “Reaction  Mechanisms  For  Toluene  Pyrolysis”  has  been  prepared 
by  M.  B.  Colket  and  D.  J.  Seery  and  has  been  accepted  for  presentation  at  the  25th 
International  Symposium  on  Combustion,  July  31  -  Aug.  5,  1994  in  Irvine  California  and 
has  been  accepted  for  publication  in  the  symposium  volume.  A  copy  of  this  paper  is 
provided  in  Appendix  A  of  this  report. 

A  paper  entitled  “Predictions  of  Soot  Dynamics  in  Opposed  Jet  Diffusion  Flames”  has 
been  written  by  R.  J.  Hall,  M.  D.  Smooke,  and  M.  B.  Colket.  It  has  been  submitted  for 
publication  in  the  special  honorific  volume  for  Professor  I.  Classman  of  Princeton 
University  as  one  of  Ae  Combustion,  Science  and  Technology  Book  Series,  published  by 
Gordon  and  Breach,  1995,  expected  publication  date.  A  copy  of  this  paper  is  in  Appendix 
B  of  this  report. 

A  manuscript  entitled  "Radiative  Dissipatitm  in  Planar  Gas-Soot  Mixtures"  authored  by 
R.  J.  Hall  and  has  been  published  in  the  Journal  of  Quantitative  Spectroscopy  and 
Radiative  Transfer,  Vol.  51,  No.  4,  pp.  635-644, 1994.  A  copy  of  this  paper  is  included 
in  this  report  as  Appendix  C. 

A  review  article  describing  state-of-the-art  modeling  of  soot  formation  processes  is  in 
preparation.  It  is  co-authored  by  I.  Kennedy  and  M.  Colket. 

rV.  Meeting  Interactions  and  Presentations 

On  Oct.  25-27,  1993,  M.  Colket  attended  the  Eastern  States  Section  of  the  Combustion 
Institute  at  Princeton  University. 

On  Oct.  28, 1993,  M.  Colket  attended  by  special  invitation  a  Colloquium  on  Combustion 
in  honor  of  Professor  Irvin  Glassman,  held  at  Princeton  University. 

On  June  1-3,  1994,  M.  Colket  attended  the  16th  Combustion  Research  Conference  held 
at  Granlibakken  in  Tahoe  City,  California  under  corporate  funds.  M.  Colket  was  invited 
to  this  DOE,  Basic  Energy  Sciences  Contractor's  Meeting  to  be  an  observer  and 
participant. 


V.  Record  Of  Inventions 


There  were  no  inventions  during  this  period. 
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Reaction  Mechanisms  for  Toluene  Pyrolysis 


ABSTRAa 

The  rich  chemistry  occurring  during  the  pyrolysis  of  toluene  has  been  investigated 
by  studying  its  decomposition  in  a  single-pulse  shock  tube  coupled  with  detailed 
chemical  kinetic  modeling  to  describe  product  formation.  This  work  provides 
information  on  the  initial  decomposition  steps  of  toluene  and  its  primary  radical  benzyl 
os  well  os  the  detailed  steps  leading  to  the  formation  of  polyaromatic  hydrocarbons. 
Temperatures  ranged  from  approximately  1200  to  1850IC  for  1%  toluene  in  argon,  total 
pressures  of  approximately  10  atmospheres  and  residence  times  near  600 
microseconds.  Pyrolysis  products  were  collected  and  analyzed  using  gas 
chrofTKJtography.  Profiles  of  hydrogen  and  hydrocarbons  ranging  from  methane  to 
pyrene  were  obtained.  This  semi-quantitotlve  data  has  been  interpreted  to  help  resolve 
many  of  the  existing  uncertainties  involving  the  pyrolytic  process.  In  particular,  this  work 
supports  arguments  for  a  low  activation  energy  process  for  benzyl  radical 
decomposition  as  evidenced  by  the  low  temperature  formation  of  cyclopentadlene 
and  the  radical  recombination  product,  benzyl-cyclopentadienyl.  Also,  the  data 
support  suggestions  for  an  important  role  of  the  methylphenyl  radical,  specifically  in  the 
formation  of  isomers  of  dimethyidiphenyl  and  as  a  key  intermediate  in  the  production 
of  anthracene.  Minimal  eviderx^e  for  rap^d  toluene  decomposition  to  phenyl  and 
methyl  was  found.  Importantly,  the  data  support  general  mechanisms  proposed  for 
ring  growth  sequences,  although  an  oddltiorKil  important  step  involving  ring  closure 
and  direct  H2  elimination  is  suggested  by  the  data.  In  addition,  the  data  indicates 
rapid  anthracene/pherranthrene  isomerizafion  above  1600  K. 
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INTRODUCTION 


The  pyrolysis  of  toluerie  and  related  chemistry  hove  been  studied  by  many 
uthors  (1-13).  The  principal  reason  for  the  more  recent  studies  has  been  that  toluene  is 
one  of  the  smallest  aromatic  molecules  implicated  in  both  soot  Inception  and  soot 
growth.  Most  of  these  studies  have  been  limited  to  the  initial  decomposition  steps  with 
little  work  on  subsequent  processes,  particularly  those  involving  polyaromatic 
hydrocarbon  (PAH)  growth  mechanisms.  Furthermore,  there  remains  a  substantial  rift 
between  European  (9,12-13)  and  US  (8.10-11)  authors  regarding  the  relative 
importance  of  the  primary  initiation  steps  In  toluene  decomposition  and  definitive 
evidence  identifying  benzyl  radical  reaction  products  is  lackir^. 

It  is  the  objective  of  this  study  to  present  some  new  data  and  interpretations  which 
will  clarify  these  primary  uncertainties  and  provide  r>ew  information  on  PAH  formation. 

Toluene  pyrolysis  studies  were  initiated  in  this  laboratory  using  a  single-pulse  shock 
tube  (SPST)  more  than  10  years  ago  (8).  Only  trace  concentrations  of  odd  carbon 
species  (except  for  methane)  were  observed,  yet  proposed  mechanisms  included 
Initial  formation  of  C5  and  C3  species  frcxn  fragmentation  of  benzyl  radicals.  In 
addition,  recovery  of  carbon  was  less  than  26%  near  temiDerotures  at  which  toluene 
decomposes  rapidly  (above  1500K).  Subsequently,  additional  work  was  performed, 
including  (1)  pyroiyzing  compounds  such  as  aiiene,  cyclopentadiene,  norbomodiene, 
and  phenylacetylene  and  (2)  collection  arKl  Identification  of  high  molecular  weight 
species  (14)  using  a  gas  chromatograph/rrx3ss  spectral  detection.  The  allene  work  was 
described  previousiy  (8)  and  helped  to  identify  rapid  formation  of  benzer^  from  C3 
species.  In  this  paper,  the  recent  results  on  the  pyrolysis  of  toluene  will  be  presented 
and  interpreted. 
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EXPERIMENTAL  APPARATUS  AND  PROCEDURES 


Experiments  were  performed  in  a  3.5  cm  singie-puise  shock  tube  using  hardware 
and  procedures  previously  described  (15).  Modifications  were  made  to  the  gas 
sampling  and  analysis  systems  to  collect  and  more  readily  identify  high  molecular 
weight  species.  The  modifications  included  higher  sample  line  temperatures  (150‘’C) 
and  use  of  a  thermal  conductivity  detector  for  hydrogen  analysis  and  a  mass  selective 
detector  for  identification.  A  flame  ionization  detector  was  used  for  quantification. 
Concentrations  of  the  larger  PAH's  reported  in  this  study  may  be  in  error  by  as  much  as 
a  factor  of  two  due  to  uncertainties  In  collection  efficiencies  and  calibration 
procedures.  More  than  50  individual  peaks  were  identified  as  products  arvj  the  more 
dominant  are  reported.  Benzyl-cyclopentadienyl  (b-cpd)  was  Identified  by  its  total 
moss  (156)  and  two  very  strong  peaks  at  91  (benzyl)  and  65  (cyclopentadienyl).  The  7 
isomers  of  dimethyidiphenyl  were  noi  irtdividually  identified.  The  data  reported  here  are 
the  combination  of  all  the  isomers  (total,  not  average).  Five  isomers  of 
diphenylmethane  at  168  amu  and  four  isomers  of  pyrene  at  202  amu  were  similarly 
combined.  The  data  presented  in  Figs,  la  to  4a  were  obtained  from  two  separate 
series  of  experiments,  although  some  of  the  data  were  averaged  to  clarify  the  figures. 

DESCRIPTION  OF  RESULTS 

Toluene  (1%  in  argon)  has  been  pyrolyzed  at  total  pressures  of  about  10 
atmospheres  for  approxirrKrtely  600  microseconds.  FirK3l  products  for  a  series  of 
individual  runs  are  presented  in  Rgs.  la-4a  for  initial  post-shock  temperatures  ranging 
from  about  1200  to  2000K.  Many  features  of  these  data  hove  been  reported  in  previous 
studies  using  a  variety  of  experimental  techniques  (1,4-6.11.16).  Methane,  acetylene, 
diocetylene.  and  benzerte  are  dominant  low  molecular  weight  products.  Ethylene  has 
been  reported  (11),  but  also  obsenred  hwe  are  vinylacetylene.  propene.  C3H4's  and 
cyciopentadier^.  Most  of  the  high  molecular  weight  species  observed  here  have 
been  reported  previously  (2.6  and  16).  Of  particular  note  is  the  new  observation  of  b- 
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cpd  presumably  formed  due  to  the  simultaneous  presence  of  benzyl  and 
cycloiDentadienyl  radicals.  It  should  be  recognized  that  some  products  (such  as  b-cpd) 
may  have  been  formed  at  least  in  part  during  the  quenching  process  in  the  SPST. 

DISCUSSION 

A  detailed  chemical  kinetic  mechanism  has  been  constructed  for  comparison 
with  these  data.  This  mechanism  is  reproduced  in  Ref  (1 7).  Because  of  space  limitations 
only  key  reactions  and  their  associated  rate  constants  are  listed  in  Table  1 .  The  reaction 
sequence  was  constructed  Initially  from  that  reported  by  Brouwer.  Muller-Markgrof,  and 
Troe  (12)  but  was  modified  substantially  to  Include  steps  describing  the  formation  of 
aromatics  and  other  species.  Also,  reactions  involving  the  rtTethylphenyl  radical  as 
earlier  suggested  by  Roo  and  Skinner  (10)  arKj  by  Pamidimukkala.  Kem.  Patel.  Wei.  arxj 
Kiefer  (11)  and  related  steps  were  itx^luded.  The  reaction  sequence  involving  light 
species  is  similar  to  those  reported  by  (16)  and  (18).  Thermodynomic  data  for  severol 
aromatic  species  were  obtained  from  the  literature  and  unavailable  thermodynamics 
were  estimated  using  group  additivity  techniques.  Thermodynamic  data  for  selected 
(aroTTKitic  and  polyaromatic)  species  are  also  provided  in  Ref  (17).  CHEMKIN  II  (20)  and 
a  modified  version  of  the  Sandia  shock  tube  code  (21)  were  used  to  perform  the  kinetic 
calculations.  Modifications  enabled  inclusion  of  experimental  pressure  decays  (along 
with  temperature  decays  assuming  isentropic.  adiabatic  expansion)  to  estimate  the 
effects  of  the  quenching  process.  Predictions  are  plotted  in  Figs.  Ib^b.  In  order  to 
understand  the  behavior  of  this  kirtetic  set.  both  a  reaction  path  analysis  and  SENKIN 
(22)  were  used. 

Identification  of  the  decomposition  pathways  of  toluene  and  its  initial  radicals 
remains  a  significant  uncertainty  in  recent  literature.  Brouwer,  et  al.  (12)  suggest  that 
k2-0.1ki,  whereas  Pamidimukkala,  et  al.  (11)  require  k2»2.5k]  in  order  to  interpret  laser 
schlleren  data.  The  primary  formation  of  ber>ryl  vs.  phenyl  leads  to  substantially 
different  decomposition  pathways.  The  conclusion  by  (11)  that  k2  is  faster  than  k^ 
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seems  to  be  unusual  in  light  of  the  endothermicities  of  the  two  reactions  (about  88  and 
102  kcal/mole  for  Reactions  1  and  2,  respectively).  The  14  kcal/mole  differential 
represents  two  orders  of  magnitude  in  the  Boitzmann  term  and  wouid  be  difficult  to 
overcome  by  differences  in  the  A-factor.  Tsang  (23)  has  recently  offered  a  iDossible 
interpretation  that  following  Reaction  1.  subsequent  thermal  equilibration  and 
recombination  of  the  products.  H+benzyl=^(C7H8)*.  leads  to  formation  of  a  thermally 

excited  toluene  molecule  with  sufficient  energy  to  decompose  to  any  of  a  variety  of 
products.  At  low  pressures,  this  process  might  be  likely  since  collisional  deactivation  is 
slow.  For  the  high  pressures  in  this  study,  the  Tsang  suggestion  appears  less  likely. 

At  low  temperatures  (-1200  to  1300K  in  this  study),  bercene.  methar>e.  hydrogen, 
and  blberucyl  are  dominant  products.  The  first  two  products  could  be  the  result  of 
Initiation  by  Reaction  2.  but  we  favor  Price's  (4)  interpretations  that  these  products  are 
due  principally  to  the  sequence  1,8,4.  whereas,  hydrogen  and  bibenzyl  are  formed 
from  Reactions  3  arKi  9,  respectively,  indeed,  reaction  pathway  and  sensitivity  analysis 
sup>port  this  interpretation,  at  least  for  the  mechanism  used  in  this  study.  In  (11). 
contributions  from  Reaction  8  were  not  Included,  since  this  reaction  is  nearly  thermally 
neutral  and  should  not  by  itself  influence  the  laser  schlieren  data. 

Use  of  the  high  value  for  k2  obtained  by  (1 1)  led  to  substantial  o  -erpredictlons  for 
methane,  ethane,  benzer^e.  and  ethylbenzene,  as  well  os  further  increosirtg  the 
otherwise  high  prediction  for  diphenyl  production  (see  Figs.  3a  and  3b).  Use  of  the 
alternative  reaction  set  (12).  with  its  high  relative  rote  for  initiation  of  toluene  to  benzyl  + 
H  leads  to  hig^  acetylene  predictions  as  shown  in  Fig.  1  b.  Much  of  this  acetylene  is  lost 
to  further  molecular  weight  growth  which  is  not  accounted  for  by  the  limited 
mechanism  of  Table  1 . 

Two  low  temperature  products  coincident  with  the  formation  of  acetylene 
(presumably  a  decomposition  product  of  benzyl  or  phenyl)  are  cyclopentadiene  (cpd) 
and  benz^yclopentodienyl  (b-cpxj)-  This  latter  product  is  almost  entirely  formed 


A-7 


during  the  quenching  process  via  the  recombination  of  the  two  stable  radicals,  benzyl 
and  cyclopentadienyl.  Traces  of  several  other  C5  species  at  masses  66  and  68  were 
observed  at  levels  typically  5  to  10  times  below  those  of  cpd.  The  presence  of  the  C5 
compounds  supp>orts  arguments  for  decomposition  of  benzyl  radical.  Cpd  is  preferred 
since  cyclopentadienyl  offers  the  lowest  energy  path  for  decomposition  (64 
kcal/mole).  The  mechanism  for  this  process  has  not  been  described  well.  Thus,  we 
explored  the  reverse  process  by  compxaring  results  from  cpd  pyrolysis  with  those  of 
cpd/C2H2  pyrolysis  (24).  The  rate  of  toluene  formation  was  enhanced  by  about  an 
order  of  magnitude  in  the  presence  of  acetylene.  Follov/ing  a  suggestion  (25),  the 
toluene  formation  may  arise  from: 

c-CsH^  +  norbomadiene  <=>  toluene  (A) 

or  c-CjHs  +  CjHj  <=>  norbomadienyl  <=>  benzyl  (B) 

Known  rates  (26)  for  (A)  do  not  account  for  the  high  temperature  dependence  (60 
kcal/mole)  of  toluene  formation  observed  in  our  cpd/C^Hj  experiments.  We  favor  (B) 
which  had  been  suggested  (27)  for  the  reverse  of  benzyl  decomposition. 
Norbomadiene  was  also  shock-heated;  only  traces  of  the  parent  remained  even  at 
HOOK  along  with  the  products,  cpd,  toluene,  cycloheptatriene,  and  acetylene, 
consistent  with  known  isomerization  and  decomposition. 

Consequently,  Reaction  10  was  adopted  as  the  preferred  low  temperature  benzyl 
decomposition  path,  A  schematic  representation  of  this  process  is  depicted  in  Fig.  5. 
The  cpd  and  b-cpd  production  were  used  to  help  determine  a  rate  for  benzyl  radical 
decomposition.  Unfortunately,  a  sensitivity  analysis  indicofes  that  the  concentrations  of 
these  species  are  dependent  on  both  the  rates  for  decomposition  of  berTzyl  (Reaction 
10)  arKj  for  cyclopentadienyl  radicals  (Reaction  13),The  rates  suggested  in  Table  1  for 
these  two  reactions  should  be  considered  accurate  to  about  a  factor  of  three  in  the 
temperature  range  of  1300  to  15CX)K.  The  activation  energies  used  here  (70  and  71 
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kcal/mole)  are  significantly  higher  than  some  values  reported  previously  (44.7  and  40, 
respectively  (12)). 

In  another  decomposition  path,  norbornadienyi  couid  isomerize  to 
cycloheptotrienyi  and  subsequently  decompose.  Alternatively,  cycloheptatriene  could 
be  formed  directly  from  a  thermally  excited  benzyl  by  bridging  of  the  CH2  group  to  the 
ortho  position,  followed  immediately  by  ring  enlargement  and  decompositton  directly 
into  two  acetylenes  and  a  propargyl  radical.  This  overall  process  requires  142  kcal/mole 
and  therefore  appears  unlikely. 

As  irKjicoted  In  Table  1,  an  alternative  mechanism  for  bercyl  decomposition  was 
included  in  part  because  of  the  results  of  (6,16  and  28)  suggesting  two  competitive 
decomposition  processes,  and  the  observations  of  vinylacetylene  in  this  study. 
Predictions  of  vinylacetylene  were  found  to  be  nearly  directly  proportional  to  the  rate 
constant  selected  for  benzyl  decomposition  to  C4H4+C3H3.  The  initial  step  for  this 
process  may  be  similar  to  that  suggested  by  Rao  and  Skinner  (29),  although  the  final 
products  suggested  here  differ  than  from  the  previous  proposal.  C4H3+C3H4  was  not 
included  since  it  is  endothermic  by  120  kcal/mole  compared  to  about  60  and  100 
kcal/mole  for  Reactions  10  and  11.  respectively.  The  rote  selected  in  Table  1  was 
selected  principally  by  the  peak  production  of  vinylacetylene. 

Other  steps  for  decomposition  of  toluene/benzyl  could  involve  the  methylphenyl 
(CH3C6H4)  radIcaKs).  Pamidumukkia.  et  al.  suggested  the  formation  of  such  radicals  by 
H-abstraction  from  the  ring.  They  then  allowed  its  subsequent  decomposition  following 
rate  kinetics  of  the  phenyl  radical.  Although  this  is  probably  a  low  limit  estimate  of  the 
decomposition  rate  constant,  the  use  of  this  rate  leads  to  a  negligible  contribution  to 
overall  ring  fracture  relative  to  the  direct  decomp)osltfon  of  benzyl  radical. 

Growth  steps  have  been  nearly  ignored  by  most  of  the  recent  kinetic  studies 
despite  substantial  documentation  of  PAH  production  at  lower  (2,3,5)  as  well  as 
comparable  temperatures  (16).  Brouwer,  et  al.  included  bibenzyl  and  trans-stilbene  in 
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their  model;  however,  the  SPST  results  indicate  the  spectrum  of  products  is  much  richer. 
Some  of  these  high  molecular  weight  species  are  products  of  combining  stable,  long- 
lived  radicals,  such  as  benzyl,  cyclopentadienyl.  and  indenyl.  partially  or  wholly  in  the 
rarefaction  wave.  Condensed  ring  aromatics  and  other  species  such  as  diphenyl  and 
diphenylmethane  (dpm)  are  principally  formed  during  the  period  of  high  temperature 
pyrolysis,  although  concentrations  can  be  altered  by  about  10%  as  the  reacting  gases 
begin  to  cool. 

Figure  3  Indicates  that  although  bibenzyl  is  the  dominant  high  molecular  weight 
species  at  low  temperature  (and  presumably  at  early  pyrolysis  tlrms).  the  combination 
of  dmdp  (dimethykjiphenyl)  isomers  and  dpm  Isomers  each  attain  nearly  comparable 
levels.  These  high  relative  levels  of  dmdp  are  comparable  to  those  obtair>ed  by  Brooks, 
et  ai.  (5).  The  dmdp  production  steps  included  in  Table  1  are  logical  extensions  of  the 
reactions  involving  the  methylphenyl  radical,  whose  possible  role  In  toluene  pyrolysis 
was  discussed  (11).  Some  dmdp  can  also  be  formed  by  benzyl  addition  to  toluene. 
Dpm  could  be  formed  by  H-atom  substitution  on  dmdp  or  benzyl  radical  addition  to 
benzene. 

As  shown  in  Fig.  4a.  several  three-ringed  species  eventually  attain  levels  similar  to 
those  of  bibenzyl,  including  fluorene,  anthracene,  phenanthrene.  and 
acenaphthylene.  Two  unusual  features  are  apparent.  Firstly,  fluorene  was  produced  at 
fairly  high  levels  for  the  low  temperature  pyrolysis.  Secondly,  anthracene  was  initially 
produced  at  levels  substantially  greater  than  its  more  stable  isomer,  phenanthrene. 
Several  ottenipts  were  made  to  model  fluorene  production.  Following  dpm  production, 
abstraction  of  an  ortho  H-atom  should  lead  to  rapid  ring  closure  and  H-atom 
ellmirK3tion.  Unfortunately,  unless  an  unreasonably  low  activation  energy  was  assigned 
to  the  H-atom  abstraction  process,  fluorene  production  was  more  than  an  order  of 
ma^ltude  too  slow.  As  an  alternative,  a  direct  four-center  elimination  of  H2  and  ring 
closure  is  proposed  (see  Fig.  6).  The  low  pre-exponential  (10®  sec’b  and  the  activation 
energy  (32  kcal/mole)  seem  acceptable.  Calculations  (30)  Indicate  that  reorientation 
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of  dpm.  so  that  the  ortho  H-atoms  are  within  ciose  proximity,  requires  oniy  5-10 
kcal/mole.  Previously.  Badger  and  Spotswood  (2)  impiicated  recombination  of  phenyi 
and  benzyl  as  a  precursor  to  fluorene  formation. 

Low  temperature  production  of  anthracene  led  to  unusually  low 
phenanthrene/anthracene  ratios:  typically  equilibrium  levels  had  been  attained  in 
other  SPST  pyrolyses  at  elevated  temperatures.  In  contrast  the  opposite  trend  was 
observed  during  phenylacetylene  pyrolysis  in  the  SPST.  In  Fig.  7.  experimental  and 
equilibrium  phenanthrene/anthracene  ratios  are  plotted  as  a  function  of  temperature. 
The  ratios  from  phenylacetylene  pyrolysis  are  near  10  at  low  temp>eratures.  but  from 
toluene  pyrolysis  are  about  0.5.  As  temp)erature  increases  above  1600K,  the  ratios  are 
both  driven  to  equilibrium  levels  near  2.5.  To  explain  the  disparate  results  at  low 
temperatures,  the  reaction  depicted  in  Fig,  6  is  proposed,  Anthracer^  is  assumed  to 
be  formed  from  a  specific  Isomer  of  dlmethyldiphenyl.  I.e.  C6H5CH2C6H4CH3  or  o- 
benzyltoluene,  where  the  attached  methyl  g^oup  is  located  In  the  ortho  position. 
Dihydroanthracene  can  be  formed  by  abstraction  of  the  benzyl  H  from  the  ortho 
methyl  group,  subsequent  ring  closure,  and  H-otom  elimination,  or  by  direct  ring  closure 
arid  H2  elimination.  Dehydrogenation  of  dIhydroanthracene  is  assumed  to  occur  at  a 
rote  comparable  to  that  for  dehydrogenation  of  1.4-cyclohexadiene  (26).  Several 
authors  have  reported  the  early  formation  of  anthracene  (3,31).  The  proposed 

mechanism  for  formation  of  anthracene  Is  similar  to  a  previous  suggestion  (32).  The 

( 

facile  conversion  of  o-benzyltoluene  has  been  demonstrated  by  Errede  and  Cassidy  (3) 
who  found  a  91%  conversion  of  o-benzyltoluene  to  anthracene. 

The  rapid  formation  of  phenanthrene  during  the  phenylacetylene  pyrolysis 
presumably  is  due  to  phenyl  addition  to  phenylacetylene  followed  by  a  cyclization 
process.  This  sequence  is  likely  similar  to  the  conversion  of  cis-stilbene  to  phenanthrene 
piLB  H2. 
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The  fastest  step  for  forming  naphthalene  is  assumed  to  be  Reaction  14 
isomerization  of  the  recombination  product  of  benzyl  and  propargyl,  There  is  little 
previous  evidence  for  this  step,  but  more  conventional  processes  (i.e..  sequential 
addition  of  acetylene  to  phenyl  (33)  and  methyl  addition  to  indene)  did  not  result  in 
reasonable  naphthalene  production  rates.  A  variety  of  steps  have  been  included  in  an 
attempt  to  describe  pyrene  formation.  By  far  the  most  important  is  phenyl  addition  to 
naphthalene  (or  naphthalenyl  addition  to  beiTzene)  followed  by  a  four-centered  H2 
elimination  and  ring  closure  to  form  fluoranthene  (see  Fig.  6).  Subsequently, 
fluoranthene  may  isomerize  to  pyrene  or  one  of  the  other  isomers  at  mass  202. 

CONCLUSIONS 

Single-pulse  shock  tube  pyrolysis  of  toluene  and  detailed  chemical  kinetic 
modeling  has  been  used  to  clarify  several  uncertain  features  of  the  decomposition 
path  and  PAH  growth  processes.  Initiation  of  the  reaction  occurs  principally  by 
decomposition  of  toluene  to  benzyl+H.  Primary  C-C  cleavage  also  occurs,  more  slowly. 
Benzyl  radical  decomposition  to  acetylene  plus  cyclopentadienyl  radical  was 
estimated  to  have  a  rate  constant  of  6x10^3  exp(-7(XXX)/RT)  sec'^  in  support  of  (12).  but 
4  times  above  the  revised  determination  (13)  at  MOOK,  Formation  of  norbomadlenyl  is 
proposed  as  an  intermediate  In  this  decomposition  process.  The  large  concentrations 
of  dimethyidiphenyl  Isomers  are  consistent  with  earlier  results  (2)  and  strongly  support 
the  suggestion  (12)  of  the  importance  of  the  methylphenyl  radical.  Mechanisms  and 
associated  rate  constants  for  forrrKatlon  of  fluorer»e  and  anthracene  are  prop>osed  and 
evidence  for  4ie  anthracene<^phenanthrene  Isomerization  is  presented.  An  estimated 
rote  constant  8x10^2  exp(-65000/RO  sec"^  is  suggested.  Ring  growth  steps  follow  similar 
paths  and  rotes  as  indicated  in  previous  studies  (33)  except  that  additional  possible 
processes  have  been  Identified.  Firstly,  four-centered,  concerted  reactions  involving 
direct  H2  elimination  and  ring  closure  with  rate  constants  of  10®  exp(-E/RT)  sec*\  where 
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30<E<40  kcal/mole.  are  suggested.  Secondly,  phenyl  addition  to  naphthalene  leading 
to  fluoranthene  formation  and  isomerization  to  pyrene  is  indicated  by  the  data. 
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TABLE  1 

PROPOSED  MECHANISM  FOR  TOLUENE  PYROLYSIS 
logk  s  logA-  E/R/T/2.303.  uNts:  cc>  mole-  sec-  kcal 

Forward  Rote 

Constant 

Reverse  Rote  Constant 

Reactions 

LogA 

E 

LogA 

E 

1  C7H8=ben2yl+H 

15.45 

88.9 

14.24 

-0.2 

2  C7H8=ph©nyl+CH3 

16.00 

97.0 

12.92 

-0.8 

3  C7H8+H=b©nzyl+H2 

14.10 

8.4 

13.50 

23.7 

4  C7H8+CH3=CH4+b©nzyl 

12.20 

11.1 

13.08 

27.1 

5  C6H6=ph©nyl+H 

16.76 

116.0 

14.46 

6.4 

6  C6H6+H=ph©nyl+H2 

14.40 

16.0 

12.71 

10.7 

7  C6H^+CH3=ph©nyl+CH4 

13.30 

16.0 

13.10 

11.5 

8  C7H8+H=C6H6+CH3 

13.08 

5.1 

12.30 

17.0 

9  2benzybbibenzyl 

12.70 

0.5 

14.93 

60.0 

1 0  b©nzyl=c-C5H5+C2H2 

13.78 

70.0 

11.57 

8.3 

1 1  benzybC4H4+C3H3 

14.30 

83.6 

8.28 

-15.6 

1 2  b©rv2yl+c-C5H5=b-cpd 

12.70 

0.0 

13.81 

49.3 

1 3  c-C5H5=C2H2+C3H3 

15.00 

71.0 

9.75 

-3.2 

14  C3H3+b©n2yl=naphthal©n©+2H 

11.78 

0.0 

0.0 

0.0 

1 6  haphthal©nyl+C6H5-4pyr©n©+H+H2 

12.00 

5.0 

0.00 

0.0 

1 6  naphthal©n©+ph©nyl-»pyr©n©+H+H2 

11.00 

5.0 

0.00 

0.0 

1 7  DPM=:fluoren©+H2 

8.00 

32.0 

10.24 

18.3 

1 8  o-b©nzyltolu©n©=dhanthrac©n©+H2 

8.00 

40.0 

9.67 

10.7 

19  anthrac©n©=ph©nanthr©n© 

12.90 

65.0 

12.73 

67.0 
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PROPOSED  PATHS  FOR  PAH  PRODUCTION 
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Data  from  phenylacetylene  pyrolysis 


A-28 


TABLE  1 

PROPOSED  MECHANISM  FOR  TOLUENE  PYROLYSIS 
logk  8  logA-  E/R/T/2.303,  units:  cc>  mole>  sec*  kcal 


Fofward  Rate 

Constant 

Reactions 

LogA 

E 

1  C7H0=benzyl+H 

15.45 

88.9 

2  C7H8=phenyl+CH3 

16.00 

97.0 

3  C7H8+H=ben2yl+H2 

14.10 

8.4 

4  C7H8+CH3=CH4+benzyl 

12.20 

11.1 

5  C6H6=phenyl+H 

16.76 

116.0 

6  C6H5+H=phenyl+H2 

14.40 

16.0 

7  C6H6+CH3=phenyl+CH4 

13.30 

16.0 

8  C7H8+H=:C6H6+CH3 

13.08 

5.1 

9  2benzyls!blbenzyl 

12.70 

0.5 

1 0  benzyls:c-C5H5+C2H2 

13.78 

70.0 

1 1  ben2yl=C4H4+C3H3 

14.30 

83.6 

1 2  ben2y1+c-C5H5=b-cpd 

12.70 

0.0 

1 3  C-C5H58C2H2+C3H3 

15.00 

71.0 

1 4  C3H3+benzyl=naphthalene+2H 

11.78 

0.0 

1 5  naphthalenyl+C6H6->pyrene+H+H2 

12.00 

5.0 

1 6  naphthalene+phenyl-4pyr©ne+H+H2 

11.00 

5.0 

1 7  DPM=fluor©n©+H2 

8.00 

32.0 

1 8  o-b©nzyltolu©ne=dhanthrac©n©+H2 

8.00 

40.0 

19  anthrac0n©=ph©nanthr©n© 

12.90 

65.0 
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Predictions  of  Soot  Dynamics  in  Opposed  Jet  DifRisicm  Flames 
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By 

R.  J.  Hall,  UTRC,  East  Hartford,  CT 
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and 

M.  B.  Coiket,  UTRC,  East  Hartford,  CT 


Abstract 


Dynamical  equations  for  spheroid  growth  have  been  coupled  to  conservation  equations  for 
opposed  jet  diffusion  flames  with  complex  chemistry  and  transport  to  provide  an  analyt¬ 
ical  tool  for  studies  of  soot  formation.  The  aerosol  dynamical  equations  use  the  discrete 
sise  spectrum,  or  sectional  method  with  inception,  surface  growth/oxidation,  and  coales¬ 
cence.  Particle  transport  by  thermophoresis  and  diffusion  is  included.  Particle  inception 
is  based  on  an  expression  involving  calculated  local  concentrations  of  acetylene  and  ben¬ 
zene  derived  from  complex  chemical  mechanisms.  Surface  growth  occurs  at  a  rate  that  is 
governed  by  temperature  and  local  concentration  of  acetylene,  and  oxidation  by  concen¬ 
trations  of  oxygen  and  hydroxyl  radical.  Thermal  radiation  by  particles  and  gas  phase 
species  is  represented  by  a  sink  term  in  the  hydrodynamic  energy  equation.  Scrubbing 
of  gas  phase  species  by  particle  growth  and  particle  thermochemistry  are  also  treated. 
The  result  is  a  strongly  coupled  gas-particulate  system;  computations  demonstrate  the 
complex  interaction  between  gaseous  species,  soot  production,  temperature  and  radiation 
and  exclusion  ot  any  of  the  above  effects  can  lead  to  significant  errors.  Comparisons  of 
model  predictions  with  experimental  data  for  methane  flames  are  presented,  and  show  ex¬ 
cellent  qualitative  agreement.  Extensive  parametric  variations  exhibiting  the  importance 
of  strain  rate,  thennal  radiation,  particulate  scrubbing  of  gaseous  species,  thermophoresis, 
and  oxidation  will  be  presented. 
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Nomenclature 


A 

At 

CptCpk 

c 

D 

D 

fv 

G 

hk 

Ib 

J 

k 

knw 

Ka,Kb,K,.KT 

m,  M 

M 

n 

N 

Na 

P 

PR 

qr 

Q 

R 

R. 

s 

t 

T 

u 

V 

VT 

VDl 

V 

Vky 

V, d,. 

Vw.» 

V 

if 

W,  W 


Strain  rate 

Particle  siirface  area  (per  particle) 

Total  particle  surface  area  per  unit  volume 

Gas  mixture  heat  capacity;  species  k  heat  capacity 

Constant  in  expression  for  soot  radiative  power  density 

Particle  diameter 

Diffusion  coefficient 

Soot  volume  fraction 

Net  surface  growth;  Sectional  surface  growth  coefficients 

Gas  species  k  enthalpy  of  formation 

Planck  function  evaluated  at  band  center 

Reduced  pressure  gradient  in  radial  direction 

Boltzmann’s  constant 

Harris- Wiener  surface  growth  rate 

Rate  parameters  in  Nagle,  Strickland-Constable  oxidation  rate  expression 

Particle  mass;  average  particle  mass  per  section 

Number  of  particle  sections 

Sectional  number  density  distribution  function 

Species  or  sectional  number  density 

Avogadro’s  number 

Partial  pressure 

Van  der  Waals  coalescence  enhancement  factor 

Net  radiative  flux 

Sectional  mass  density 

Speciflc  surface  growth  or  oxidation  rate 

Gas  constant 

Mass  source/sink  rate 

Time 

Gas  temperature 
Raffial  flow  velocity 
Axial  flow  velocity 
Thermophoretic  velocity  of  sections 
Diffusive  velocity  of  section  1 
Axial  mass  flux 

Diffusive  velocity  of  gas  species  k 
Fuel  stream  mass  flux 
Oxidize  stream  mass  flux 
Thermal  collisional  velocity  of  sections 
Molar  production  rate 

Species  molecular  weight;  average  molecular  weight 
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H  ^ 


r  Radial  coordinate 

Axial  coordinate 

Mass  fraction  of  species  or  section 

Greek  symbols 

a  Integrated  band  intensity;  thermal  accommodation  coefficient,  Frenklach-Wang 

surface  growth  steric  factor 

Coalescence  collision  frequency;  inter-sectional  average  collisional  frequency 
S  Dirac  delta  function 

A  Thermal  conductivity 

^  Gas  viscosity 

u  Stoichiometric  coefficient  in  gas-solid  scavenging  analysis 

(  Mixture  fraction 

PtPtiPi  Gas  density;  soot  specific  density;  density  of  i-th  gasesous  species 

x'  Parameter  in  Nagle,  Strickland-Constable  oxidation  rate  expression 


I.  Introduction 

A.  Effects  of  soot  formation  in  flames 

Many  of  the  effects  of  soot  formation  have  been  discussed  in  detail  in  a  great  number  of 
technical  and  review  papers  (see,  e.g.  Wagner,  1979;  Haynes  and  Wagner,  1981;  Glassman, 
1988;  Smith,  1981;  and  Barfknecht,  1983).  These  articles  have  discussed  pollution  (includ¬ 
ing  carcinogens)  plume  visibility,  component  lifetimes,  and  radiation.  More  recently,  as 
combustion  models  have  advanced,  the  significant  importance  of  soot  formation  to  flame 
modeling  has  also  been  recognized. 

The  two  most  obvious  and  acknowledged  flame  processes  affected  by  soot  are  radiation 
and  emissions.  Fractional  energy  losses  from  radiation  as  high  as  twenty  to  thirty  percent 
(Markstein,  1984;  Hall  and  Bonczyk,  1991)  are  not  unusual  for  coflow  diffusion  flames. 
Emissions  have  long  been  known  to  effect  visibility  of  combustion  exhaust  and,  most  re¬ 
cently,  signicant  health  effects  of  combustion-generated  soot  have  been  highlighted  in  a 
New  York  Times  front  page  article  (July  19,  1993).  The  depression  of  temperature  due 
to  radiative  losses  can  in  turn  effect  flame  length  (Bhattacharjee  and  Grosshandier,  1988) 
and  other  temperature-dependent  processes  such  as  formation  of  other  pollutants,  partic¬ 
ularly  NOx.  A  secondary  but  potentially  important  effect  in  a  strongly  sooting  flame  is 
the  fuel/soot  thermochemistry  which  can  vary  dramatically  depending  on  fuel-type.  Also, 
soot  particles  formed  in  a  pyrolytic  region  can  then  be  transported  into  an  oxidizing  zone 
where  oxidation  processes  are  much  slower  than  the  kinetics  of  gas-phase  oxidation  thus 
altering  the  heat  release  rates  and  the  temperature  profiles.  In  flames  in  which  substantial 
fractions  of  fuel  carbon  are  converted  to  soot,  this  local  extraction  of  fuel  carbon  can  lead 
to  an  alteration  of  the  effective  local  equivalence  ratio  of  the  gas-pahse  species  leading  to  a 
shift  in  the  local  H2/H3O  and  GO/CO3  conversion  fractions  and  hence  local  heat  release. 
Thus,  for  flames  in  which  flame-front  chemistry  controb  pollutant  production  (e.g.  prompt 
NO)  these  many  temperature-perturbing  effects  of  soot  formation  may  be  non-trivial. 


B.  Brirf  review  of  soot  formation  modeling 

During  the  past  several  years,  there  has  been  a  dramatic  increase  in  attempts  at  modeling 
soot  formation  in  premixed  and  both  laminar  and  turbulent  diffusion  flames.  By  neces¬ 
sity,  however,  as  flame  geometry  and  structure  increases  in  complexity  the  assumptions 
r^arding  the  soot  model  become  more  simplified.  So  far,  reasonable  agreements  have  been 
obtained  for  many  of  the  flames  examined.  These  excellent  results  should  be  considered  a 
major  success  story  in  combustion  modeling.  For  laminar,  premixed  flames,  Frenklach  and 
coworkers  (1985,  1987, 1988  and  1990)  are  using  the  most  sophisticated  soot  model  which 
includes  detailed  reaction  chemistry  up  to  and  including  large  PAH  structures,  perhaps 
the  most  realistic  inception  concepts,  detailed  growth  and  oxidation  steps,  and  conven¬ 
tional  aerosol  djmamics  including  size-dependent  processes.  This  approach  has  been  very 
successful!  at  predicting  sooting  limits  in  premixed  flames  (Markatou,  et  al,  1993).  Colket 
and  Hall  (1994)  have  used  generally  similar  approaches  for  premixed  flames  although  a 
sectional  model  has  been  used  to  simplify  the  aerosol  equations.  In  addition,  they  related 
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the  inception  species  current  to  the  benzene  current  using  the  method  of  characteristics. 
Although  this  assumption  overestimates  inception  and  is  physically  unrealistic,  it  allows 
for  a  great  simplification  in  the  gas-phase  kinetics  while  still  treating  the  general  features 
of  the  stability  of  aromatics  and  its  importance  to  soot  production. 

For  more  complex  fiames  greater  simplifications  have  been  utilized,  although  there  is  some 
concern  for  extrapolation  of  these  models  to  significantly  different  conditions,  i.e  use  of 
practical  fuels,  high  inlet  temperatures,  high  pressures,  etc.  These  more  global  models  use 
simpler  inception  concepts  and  assume  a  monodisperse  size  distribution  with  two  dynam¬ 
ical  equations,  one  for  the  particle  number  density  and  the  other  for  the  volume  fraction. 
Specific  concerns  include  the  use  by  Kennedy,  et  al.  (1990)  of  empirical  specific  soot  growth 
rates  (expressed  as  a  function  of  mixture  fractions)  obtained  for  the  same  fuel  (ethylene)  in 
an  opposed  jet  diffusion  flame  operating  at  the  same  ambient  conditions  as  for  their  model 
of  coflow  laminar  and  turbulent  diffusion  flames.  Predictions  were  found  to  be  essentially 
independent  of  assumed  inception  rates.  Others  assume  inception  is  proportional  to  the 
acetylene  concentration  and  a  surface  growth  rate  proportional  to  the  square  root  of  the 
surface  area.  Sivathanu  and  Gore  (1993)  have  shown  this  approach  provides  good  results . 
for  coflow  diffusion  flames  while  Fairweather,  et  al.  (1992)  and  Leung,  et  al.  (1991)  have 
incorporated  these  assumptions  into  a  flamelet/p.d.f.  code  for  turbulent  jets.  Extrapola¬ 
tion  of  these  assumptions  to  ambient  conditions  far  outside  the  range  of  tested  conditions 
is  probably  uncertain.  Typically  for  these  codes,  radiative  loss  fractions  are  assumed.  As 
Sivathanu  and  Gore  have  shown,  however,  there  can  be  a  strong  coupling  to  radiation. 
In  their  study,  changes  in  the  radiative-loss  fraction  of  just  15%  were  found  to  effect  dra¬ 
matically  temperatures  and  hence  soot  kinetics.  Consequently,  predictions  of  soot  volume 
fractions  were  altered  by  as  much  as  an  order  of  mi^itude. 


C.  Objectives  of  this  work 

The  objectives  of  this  work  are  to  develop  a  code  which  fully  integrates  detailed  gas-phase 
chemical  kinetics  with  known  processes  involved  with  soot  pr'^duction  in  flames,  includ¬ 
ing  inception,  growth,  oxidation,  agglomeration,  radiation,  thermochemistry,  scavenging, 
and  transport.  We  have  selected  a  laminar,  opposed  jet  flame  code  (Smooke,  Puri,  and 
Seshadri,  1986)  for  modification.  An  advantage  of  such  a  configuration  is  that,  once  the 
gas-phase  chemistry  is  defined,  solutions  could  be  obtained  for  a  variety  of  flames  of  dif¬ 
fering  fuel  type  and  ambient  conditions  from  which  ‘empirical’  growth  rates,  etc.  could 
be  determined  as  input  to  codes  describing  more  complex  flame  geometries.  It  is  not  the 
purpose  of  this  study  to  identify  specific  or  recommended  flame  chemistry  for  formation 
of  aromatics  (or  PAHs),  or  for  selecting  valid  inception  or  growth  models.  Rather,  it  is 
an  objective  of  this  study  to  develop  a  tool  (computer  code)  which  is  general  enough  such 
that  a  variety  of  such  models /concepts  could  be  examined  and,  consequently,  sensitivities 
to  uncertainties  of  such  processes  can  be  estimated. 

This  manuscript  describes  details  of  the  code  and  assumptions  used  in  its  development. 
In  addition,  perturbation  effects  due  to  inclusion/exclusion  of  varioxis  phenomena  such  as 
oxidation,  species  scavenging  (by  soot)  and  radiation  from  soot  or  from  gas-phase  species 
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are  examined.  Methane  is  assumed  to  be  the  fuel  for  these  studies.  Since  methane  is 
a  relatively  lightly  sooting  fuel,  many  of  the  effects  described  herein  are  probably  low 
limit  estimates.  As  part  of  this  effort  a  kinetic  analysis  compeu’ing  various  routes  to  the 
formation  of  benzene  in  a  methane  fired  system  has  been  performed  in  order  to  minimize 
additional  reactions  and  species  required  to  model  benzene  formation  and  destruction.  The 
major  conclusions  of  this  investigation  are  included. 


II.  Description  of  Model 

A.  Problem  formulation  -  opposed  jet  diffusion  flame  code 


To  model  the  combined  gas  phase/soot-laden  counterfiow  diffusion  flame  (Figure  1),  we 
utilize  the  elliptic  form  of  the  two-ulmensional  conservation  equations  in  cylindrical  coor¬ 
dinates  for  the  gas  phase  together  with  sectional  species  equations.  The  complexity  of  the 
entire  system  is  reduced  (Smooke  and  Giovangigli,  1992)  by  seeking  a  similarity  solution  of 
the  form  u  =  rU(z),  v  =  v(z),  Yk  =  Yk(z),  k  =  1,2,. . .K,K -h  1, . . .  ,K -f- M,  T  =  T(z), 
where  x  and  y  denote  the  independent  radial  and  axial  spatial  coordinates,  respectively;  T, 
the  temperature;  Yk,  the  mass  fraction  of  the  k^**  gas  phase  species  or  the  (k  -  K)**'  soot 
size  class  and  u  and  v  the  radial  and  axial  components  of  the  velocity.  If  we  substitute  the 
expressions  for  the  two  velocities  into  the  continuity  and  the  momentum  equations,  we  find 
that  the  reduced  pressure  gradient  in  the  radial  direction  is  constant,  i.e.,  r  ^  where 
P  denotes  the  pressure.  If  we  use  this  result,  together  with  the  similarity  assumptions, 
we  find  that  the  two-dimensional  set  of  governing  equations  can  be  reduced  to  a  nonlin¬ 
ear  two-point  boimdary  value  problem  in  the  axial  direction.  The  governing  equations  for 
mass,  momentum,  chemical  species,  soot  size  classes  and  energy  cem  be  written  in  the  form 


dz 


-H  2pU  =  0 


(1) 


and 


d  /  dU\  ,  dU  ^ 

d  dYk 

-^(pYk(vT  +  VDj)-V-^-HQk  =  0,  k  =  K-Hl,K-h2,...,K-l-M 

S  ("f )  -  "'"f  - 1  -  tw  +  *l)W.h.  +  ^  =  0 

'  '  k=l  k=l 


PW 

RgT 


(2) 

(3) 

(4) 

(5) 

(6) 


where  V  =  pv,  and,  in  addition  to  the  variables  already  defined,  p  denotes  the  mass  density; 
Wk,  the  molecular  weight  of  the  k^^  species;  W,  the  mean  molecular  weight  of  the  mixture; 
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Rg,  the  universal  gas  constant;  A,  the  thermal  conductivity  of  the  mixture;  Cp,  the  constant 
pressure  heat  capacity  of  the  mixture;  Cpk,  the  constant  pressure  heat  capacity  of  the 
species;  w{|,  the  molar  rate  of  production  of  the  species  per  unit  volume  due  to  gas 
phase  chemistry,  w}^,  the  molar  rate  of  production  of  the  k^**  species  per  unit  volume  due 
to  scrubbing/replenishment  by  the  soot  growth/oxidation  processes;  Qk  the  rate  of  change 
of  section  k  due  to  aerosol  processes;  hk,  the  specific  enthalpy  of  the  k^^  species;  /x  the 
viscosity  of  the  mixture;  Vky,  the  diffusion  velocity  of  the  k***  species  in  the  y  direction; 
vt  and  vd^  the  thermal  and  mass  diffusion  velocities  of  the  k'^^  soot  size  class  and  ^  the 
divergence  of  the  the  net  radiative  flux  which  in  the  optically  thin  limit  is  given  by  (Hall, 
1993  and  1994) 

=  CfvT®  +4jr^  Oik  Pkibik  (7) 

^  ik 

where  the  ik  summation  is  taken  over  the  optically  active  bands  i  of  the  k^*^  gas  phase 
species  and  where  C  is  a  constant  that  is  a  function  of  the  soot  index  of  refraction.  At 
low  strain  rates,  ~5-20%  conversion  of  flame  enthalpy  release  to  radiative  loss  can  occur 
(Vranos  and  Hall,  1993;  Hall,  1994),  leading  to  gas  temperature  reductions  large  enough  to 
affect  soot  formation,  as  will  be  seen.  Soot  contributions  to  the  density  and  heat  capacity 
are  usually  small  and  have  been  ignored.  Thermal  equilibrium  between  the  gas  and  the 
particulate  phases  has  also  been  assumed.  The  form  of  the  chemical  production  rates  and 
the  gas  phase  diffusion  velocities  can  be  found  in  detail  in  Giovangigli  and  Darabiha  (1988). 

The  boundary  conditions  and  solution  algorithm  for  the  resulting  two  point,  boundary 
value  problem  posed  by  Equations  1-7  are  discussed  in  Appendix  A. 

B.  Gas-phase  chemistry 

Except  for  the  aromatic  fomation/destruction  chemistry,  the  methane  reaction  set  used  in 
this  study  is  identical  to  that  used  previously  (Smooke,  et  al.,  1992)  to  examine  methane 
cofiow  flames.  This  reaction  mechanism  includes  a  relatively  small  C2-reaction  chemistry 
set.  Several  reactions  and  species  have  been  added  to  this  mechanism  to  include  formation 
and  destruction  of  benzene.  The  added  reactions  are  listed  in  Table  I.  The  entire  mecha¬ 
nism  contains  33  species  and  102  elementary  reactions.  The  benzene  kinetics  were  based 
on  an  evaluation  of  benzene  formation  mechanicms  (Colket  and  Smooke,  1993)  during  ox¬ 
idative  pyrolysis  of  methane  in  a  single-pulse  shock  and  a  recent  proposal  (Bittker,  1991) 
for  the  oxidation  of  benzene.  In  the  shock  tube  study,  six  different  steps  forming  benzene 
were  compared  including  several  C4  +  C3  reactions.  By  far,  the  most  dominant  step  was 
the  propargyl-propargyl  recombination  recently  suggested  to  be  the  predominant  benzene 
forming  step  during  allene  pyrolysis  (Kem,  et  al.,  1988)  and  in  acetylene  flames  (Miller 
and  Melius,  1992).  Modeling  using  this  reaction  predicted  benzene  profiles  in  reasonable 
agreement  with  the  single-pulse  shock  tube  data  on  benzene  production  during  partial 
oxidation  of  methane. 


TABLE  I 

SET  OF  REACTIONS  AND  RATE  COEFFICIENTS 
DESCRIBING  BENZENE  FORMATION  AND  DESTRUCTION 
logk  =  logA  +  nlogT  -  E/R/T/2.303  * 


REACTIONS 

FORWARD 

RATE 

CONSTANT 

REVERSE 

RATE 

CONSTANT 

1 

C3H4+H=C2H2+CH3 

LogA 

13.60 

2.4 

LogA 

12.28 

12.6 

2 

C3H4+H=C8H3+H2 

12.00 

1.5 

11.02 

22.4 

3 

CH3+C3H4=C8H3+CH4 

12.30 

7.7 

12.83 

29.5 

4 

2C3H3=>C6H8+H 

13.00 

0.0 

0.00 

0.0 

5 

CeHe+H^CeHs+Hj 

14.40 

16.0 

12.67 

10.5 

6 

CeHs + H=CeH6 

13.70 

0.0 

16.03 

109.9 

7 

C6H8+CH4=C6H6+CH3 

12.48 

5.0 

12.69 

9.6 

8 

CeHs^l-CeHs 

13.65 

72.5 

10.02 

8.6 

d 

1-C6H8=>2C2H2+C2H 

14.00 

36.0 

0.00 

0.0 

10 

C«H«+02*>c-C8H6+C0+0 

12.32 

7.5 

0.00 

0.0 

11 

c-CsHe^H+c-CsHg 

15.30 

81.0 

14.32 

1.2 

12 

H+C-C8Hs~H3+C-C8H8 

12.48 

8.0 

12.10 

32.5 

13 

H+c-C8He=C3H4+C8H8 

12.70 

18.0 

8.62 

14.1 

14 

CH8+c-C8Hs=CH4+c-C8H8 

12.70 

5,0 

13.83 

30.4 

15 

C-C8H8  =C2H2  +C8H8 

14.00 

74.0 

8.69 

3.6 

16 

C6He+OH=C«H8+H20 

13.33 

4.6 

12.27 

14.4 

17 

C6H6+0=C-C8H8+C0+H 

13.44 

4.9 

11.23 

1.1 

18 

C-C8H8+0=>C2H8+C2H3+C0 

13.70 

0.0 

0.00 

0.0 

19 

C6H8+OH=>c-C8H8+CO+H 

13.00 

0.0 

0.00 

0.0 

*  NOTES: 

Units  for  A:  cc,  moles,  sec. 

Units  for  E:  kcal/mole. 

«  _  *»  represents  forward  and  reverse  directions  included  in  model. 

*  -  ”  represents  forward  direction  only  included  in  model.. 
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With  this  mechanism,  profiles  of  gas-phase  species  were  calculated  using  the  conventional 
version  of  the  opposed-jet  flame  code.  Calculated  values  are  shown  in  Figures  2  and  3  for 
two  different  strain  rates.  Species  are  plotted  versus  a  the  normalized  coodinate  used  by 
Zhang,  et  al.  (1992): 


f  Z  —  Zstag 
Z  =  - — 


(8) 


where  yn  is  the  location  of  the  peak  fiame  temperature  and  yst»g  Is  the  location  of  the 
stagnation  plane.  These  predictions  can  be  compared  with  those  obtained  by  Zhang,  et 
al.  The  benzene  profile  (both  shape  and  magnitude)  as  well  as  most  other  species  profiles 
are  predicted  well  by  the  kinetics  model,  although  acetylene  concentrations  appear  to  be 
overpredicted.  Since  the  local  temperature  for  benzene  formation  is  about  1400-1500K 
in  the  flame,  while  benzene  production  in  the  shock  tube  occurred  at  temperatures  of 
1700-1800K,  it  is  concievable  that  an  alternative  benzene  formation  mechanism  may  be 
contributing  in  the  flame  system.  Consequently,  a  test  was  performed  by  including  all 
benzene  steps  (see  Colket  and  Smooke,  1993)  and  appropriate  steps  for  the  formation  and 
destruction  of  the  intermediate  species.  Using  this  larger  reaction  set,  virtually  no  change 
in  the  benzene  production  was  observed,  while  if  the  propargyl  recombination  step  was 
removed,  peak  benzene  profiles  were  reduced  by  more  than  two  orders  of  magnitude.  The 
importance  of  this  reaction  to  benzene  production  in  a  methane  system  has  also  been  found 
by  Lindstedt  and  Skevis  (1994). 

Besides  the  reduced  oxidation  sequence  for  decomposing  benzene,  an  important  step  in 
defining  the  benzene  profile  is  the  thermal  decomposition  of  the  phenyl  radical.  Changes 
in  this  rate  constant  led  to  nearly  comparable  changes  in  the  peak  concentration  of  benzene. 


C.  Soot  Inception 

A  variety  of  assumptions  for  soot  inception  have  been  made  in  the  literature.  These  include 
inception  based  on  acetylene  concentrations,  benzene  production,  as  well  as  production 
of  polyaromatic  structures  of  different  sizes.  In  order  to  minimize  complications  to  the 
reaction  chemistry  while  attempting  to  approach  the  use  of  a  more  realistic  inception 
process,  simplified  steady-state  assumptions  were  used  to  estimate  the  formation  of  two- 
and  three-ringed  aromatic  species.  Consider,  for  example,  the  two  sequences: 


CjHj  +  CsHs  CgHe  +  H 
H  +  CgHs  GgHg  +  Hj 
C3H3  +  GgHs  C10H7 


(20)  K«,,3o  = 

121}  K,q.3i  =  10^«®e-^«®VT 

(22)  k33  =  cc/mole/sec 


(9) 


or 

2C3H3  +  CsHs  -♦  CioH?  +  H2  (overall) 
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and 


CeHe  +  CoHs  ^  C13H10  +  H  (23)  Keq.23  = 

H  +  C12H10  C12H9  +  H2  [24)  Keq,24  =  jQl-36g— 2768K  (10) 

C2H2  +  C12H9  — ♦  G14H10  +  H  (25)  k25  =  cc/mole/sec 

or 

C2H2  +  CeHo  +  CeHs  — ♦  Ci<iHio  +  H  +  H2  (overall) 

Assuming  steady  state  values  of  CgHs,  GgHe,  Gi2Hio,  and  C12H9  and  that  H2  <<  C2H2, 
the  rate  of  production  of  the  polyaromatic  species  can  be  estimated  to  be 


dlCioHr) 

dt 


Keq,2oKeq,2lheq,22 


[C2H2)" 

|H»| 


|CeH,| 


and 


=  K.,,2sK«,,Mk«,,j5l^l|C«H,||C.H5l 


(12) 


or  using  the  rate  and  equilibrium  ‘constants’,  these  expressions  become 


diCioHr)  _  --j11.88  (-4S78/T)  1^2^21^  tp  TT  1 

'  iH^r' ' 


(13) 


and 


5I2^S£l  =  10‘>«'.‘-“*>/'"'l^lC.He||C,H,|  (14) 

where  the  gas  phase  concentrations  and  temperatures  can  be  evaluated  at  local  conditions. 
With  the  further  assumptions  that  inception  b  limited  by  the  formation  of  poiyaromatics 
and  Qxidation/decomposition  of  such  species  can  be  neglected,  the  inception  rate.  Si,  in 
grams/cc/sec  was  initially  assigned  to 


Si  =  127  X 


d[CioHy] 

dt 


+  178  X 


d[Ci4Hio) 

dt 


(15) 


where  the  constants  (molecular  weights)  are  provided  to  convert  from  molar  to  mass  units. 
As  will  be  indicated  later,  soot  production  using  this  rate  was  small  and  the  inception  rate 


B-11 


for  most  of  this  study  was  assumed  to  be  eight  times  that  calculated  from  the  sum  of  the 
above  terms.  For  all  cases  examined  in  this  study,  surface  growth  (Equation  lb),  was  the 
dominant  term. 

The  ‘requirement’  for  an  augmented  inception  rate  could  be  due  in  part  to  several  other 
model  uncertainties.  For  example,  the  rate  of  naphthalene  and  phenanthrene  formation  as 
calculated  by  Eqn.  (15)  could  itself  be  low.  The  recent  results  by  Colket  and  Seery  (1994) 
and  Lindstedt  and  Skevis  (1994)  suggesting  that  other  steps  contribute  significantly  to  PAH 
growth  support  this  argument.  Alternatively,  the  rate  of  surface  growth  could  be  in  error. 
Kennedy  and  coworkers  had  to  utilize  an  empirical  surface  growth  rate  term  since  the  HW 
expression  inadequately  described  the  rate  of  soot  growth  in  their  flames.  A  high  inception 
rate  could  be  required  to  offset  a  low  surface  growth  rate.  Also,  inaccuracies  in  the  kinetics 
set  could  be  in  part  responsible.  Either  the  benzene  concentration  is  underpredicted,  or 
the  oxygen  concentration  is  overpredicted  by  the  kinetics  model.  It  seems  unlikely  that  the 
benzene  concentration  is  that  poorly  predicted.  The  substantial  non-zero  concentration  of 
oxygen  on  the  fuel  side  (see  Figures  2  and  3)  seems  unusual,  although  comparable  behavior 
of  oxygen  was  observed  using  more  complete  methane  oxidation  mechanism  (Glarborg,  et 
al.,  1986)  as  well.  With  a  lower  level  of  oxygen  in  the  soot  growth  region  or  a  shifting 
of  the  oxygen  profile  away  from  the  fuel  rich  region,  predicted  soot  concentrations  can  be 
substantially  increased  without  augmenting  the  inception  rate. 

Because  of  the  potential  uncertainties  in  the  model,  the  reader  is  reminded  that  the  focus 
of  this  paper  is  the  description  of  the  model  and  the  demonstration  that  soot  formation 
is  a  result  of  complex  interactions  between  chemical  and  physical  phenomena  causing 
perturbations  to  local  temperatures  and  gas-phase  species. 


D.  Soot  Growth  Model(s) 

Several  procedures  for  calculating  the  rate  of  surface  growth  of  soot  particles  were  con¬ 
sidered.  These  include  the  Harrb  and  Weiner  (1983)  expression  (HW),  the  Frenklach  and 
Wang  (1990)  procedure  (FW)  and  a  modification  (MODFW)  of  the  FW  method  (Colket 
and  Hall,  1994).  In  this  study,' the  multiplicative  parameter,  a,  required  for  the  FW  pro¬ 
cedure  is  assigned  a  value  of  1.0,  although  this  value  may  be  an  order  of  magnitude  too 
high  near  the  high  temperatures  in  a  diffusion  flame.  For  most  of  the  results  presented 
here,  the  HW  expression 


^  =  kHW  A  Pc,H,  (16) 

has  been  used.  A  is  the  surface  area  of  the  soot  particle  and  PcaH,  b  the  partial  pressure  of 
acetylene  in  atmospheres.  Using  an  overall  activation  energy  for  the  soot  growth  process 
of  31.8  kcal/mole  (Hura  and  Glassman,  1988)  and  Harrb  and  Weiner’s  measured  rate 
at  1650K  of  3  xld®  grams/sec/cm* /atih,  knw  b  determined  to  be  47  exp(-16004/T) 
gram/sec/cm^/atm.  Ageing  (Dash,  1985)  which  reduces  particle  surface  growth  rates 
with  increasing  time  was  not  included  in  thb  problem. 
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E.  Soot  Dynamics 
1.  Aerosol  Model 

The  growth  of  soot  spheroids  has  been  modeled  as  a  classical  aerosol  dynamics  problem, 
involving  the  division  of  the  size  range  of  interest  into  discrete  intervals  or  classes  with  a 
logarithmic  transformation  of  the  size  range,  and  then  deriving  a  source  for  the  size  class 
mass  densities  with  terms  representing  inception,  surface  growth  (or  oxidation),  and  coag¬ 
ulation  (coalescence).  The  sectional  analysis  is  dbcussed  by  Gelbard  and  Seinfeld  (1980) 
and  Gelbard,  Tambour,  and  Seinfeld  (1980),  and  uses  certain  algorithms  of  the  well-known 
MAEROS  program  (Gelbard,  1982).  It  is  assumed  that  surface  growth  and  coalescence 
have  the  free  molecule  form,  and  that  the  spheroids  consist  of  the  single  component,  carbon. 

In  the  sectional  analysis,  it  is  assumed  that  the  boundaries  of  the  sections  vary  linearly 
on  a  log  scale.  Thus,  if  we  have  M  sections,  with  diameter  minima  and  maxima  D(0)  and 
D(M),  respectively,  the  sectional  boundaries  will  be  given  by 


Dt  =  D(0) 


«/M) 


(17) 


Section  t  (1<  ^  <  M)  will  thus  have  the  boundaries  m(^-l)  and  m(£)  in  particle  mass  where 


irD® 

m-*.. — 


(18) 


and  the  soot  density  p,  is  taken  to  be  that  of  solid  carbon,  1.8  g/cc.  If  n£(m,t)  represents 
the  number  density  per  u:  mass  of  particles  within  a  section,  the  total  mass  within 
the  section  is  given  by  the  integral 


m  n«(m,t)  dm 


Choosing,  as  in  Gelbard,  the  form 


(19) 


ii((m,t) 


(20) 


where 


At  =  (tn  —)  '  =  (X,  -  Xt-,)-' 

\  mr_i  / 

X  *  to(m) 
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results  in  the  following  expression  for  the  total  number  density  of  particles  in  class  i 


N«=  /“'  nt(m,t)  dm=  Q<  A,  - !-]  = 

\^e-i  rntj  me 

where  the  average  particle  mass  in  section  I  is  given  by 


/Xl,  nt(m,t)dm 


(22) 


(23) 


The  per  particle  rate  of  growth  due  to  surface  mass  addition  and  oxidation  has  been 
assumed,  as  discussed,  to  be  proportional  to  particle  surface  area  (free  molecule  or  large 
Knudsen  number  form) 


^  =  G'(t)A  =  G(t)mVs 


(24) 


where  A  is  the  particle  surface  area,  and  the  specific  growth  rate  has  the  overall  form 


G'(t)  =  (Growth  rate  by  acetylene  or  other  growth  species  addition 
—  oxidation  rate  by  O2,  OH) /unit  surface  area. 


(25) 


or'  G^(t)  =  Rg  ~  Roa  ~  Roh 

The  mass  addition  term  Rg  (see  Eqn.  16)  has  been  derived  for  several  kinetic  models 
of  the  surface  growth  process,  as  discussed  in  the  preceding  section.  Oxidation  of  soot 
by  OH  radicals  is  assumed  to  proceed  at  a  gas  kinetic  coliuion  frequency  multiplied  by 
a  collision  probability  of  0.13  (Neoh,  Howard,  and  Sarofim,  1981).  Thus,  with  Noh>  and 
Na  representing  the  OH  number  density  and  Avogadro*s  number,  respectively,  the  OH 
oxidation  is 


Roh  =  (0.13)  x  NoHy 


RgT  12 

— - —  X  — 


27rWoH  Na 


=  16.7 


POH 

Vf 


(26) 
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where  Poh  the  OH  partial  pressure  in  atmospheres,  and  the  specific  growth  rate  is  in 
c.g.s.  units.  For  oxidation  by  O3,  the  Nagle  &  Strickland-Constable  (1963)  expression  is 
used.  For  an  O3  partial  pressure  Po,  in  atmospheres,  the  specific  oxidation  rate  is  thus 

•  (1+‘k°p1)  +  KbPo,  (1  -  x'))  (27) 

where  K*,  Kb,  Kt,  Kt*  &nd  x'  as  given  by  Nagle  &  Strickland-Constable. 

At  atmospheric  pressure,  the  coalescence  rates,  0  for  flames  are  in  the  free-molecule  regime, 
and  thus  given  by 


=  PR  (m,m07(Dm  +  Dm')*  Vmm' 

4 

with 


V 


mm'  — 


(28) 


and  provision  is  made  for  a  Van-der  Waals  enhancement  factor  PR.  For  collision  partners 
of  low  and  similar  mass,  theory  predicts  enhancements  of  the  gas  kinetic  rate  by  as  much 
as  2.2  (Harris  and  Kennedy,  1988);  we  use  a  nominal  value  of  1.5  for  all  interactions.  In 
principle,  it  should  be  possible  to  include  the  mass  dependence  of  the  PR  in  the  sectional 
0  calculation. 

The  dynamic  balance  equation  for  the  can  be  represented  as 


dt 


t-i 


5  E  E  QiQi-Q«E 


1=1  j=l  i=l 

5  - «!«  E 


+  Q< 

+  *Gr+/H  G«  Qr  +  ^«Si(t) 


==  -1  G(t)  >  0 

=  -HI  G(t)  <  0 


(29) 


where  the  0  represent  sectional  coagulation  coefficients,  are  intra-sectional  growth 
coefficients,  ^G^  are  inter-sectional  growth  coefificients,  and  Si  (see  Eqn.  15)  is  a  particle 
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inception  rate  (mass  units)  for  the  initial  size  class.  The  value  of  the  index  In  is  determined 
by  whether  the  net  growth  is  positive  or  negative,  and  it  is  understood  that  the  term 
containing  it  will  not  be  present  for  ^=1  for  net  growth,  and  will  not  be  present  for  ^=M 
for  net  oxidation  (G(t)  greater  or  less  than  zero,  respectively).  Implicit  is  the  assumption 
of  complete  particle  coalescence  after  collision,  with  no  aggregate  or  chain  formation. 

Because  in  the  free  molecule  regime,  both  the  surface  growth  and  coagulation  rates  are 
factorable  into  products  of  functions  of  mass  times  functions  of  time,  a  considerable  sim¬ 
plification  of  the  task  of  calculating  sectional  coefficients  is  made  possible  because  the 
mas8>dependent  calculations  need  be  performed  only  once,  with  the  time  dependence  im¬ 
posed  during  the  course  of  the  integration  (e.g.  v/T).  With  the  general  functional 

form  for  the  surface  growth,  the  intra-sectional  growth  coefficients  can  be  calculated  ana¬ 
lytically  as 


■G«  =  3A«G{t)(^-^)  (30) 

Vm4i  m/  / 

and  the  inter-sectional  coefficients  are  related  to  these  by 

=  ^G^  -  G(t)  >  0 

m<+i  -  mi 

^Gi  =  ^Gi  G(t)  >  0 

m^-i  -  m< 

The  inter-sectional  coefficients  are  derived  from  a  number-  and  mass-conserving  algorithm 
due  to  Warren  and  Seinfeld  (1985).  If  some  other  dependence  of  surface  growth  on  A  is 
assumed  (e.g.  «  v^;  Furweather,  et  al,  1992)  a  different  but  easily  derived  expression 
for  the  ^Gi  would  be  needed.  The  sectional  coalescence  coefficients  {P)  are  calculated  from 
the  firee-molecule  collision  frequencies,  (/9)  using  standard  numerical  quadratures. 

With  the  solution  for  the  Qr,  a  number  of  aerosol  properties  can  be  evaluated.  The  soot 
volume  fraction  will  be  given  by 


(31) 


and  its  rate  of  change  by 


U  =1.  T  ^GiQi  =  -G'(t)AT 

p.  V  p> 


(32) 


(33) 
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where  At  is  the  total  particle  surface  area  per  unit  volume  (defined  below).  The  moments 
of  the  diameter  distribution  are  expressible  as 


<  D"  >= 


/  D°(m)nf(m,t)dm 
/  n<(m,t)dm 


(34) 


For  comparisons  with  particle  diameters  obtained  from  laser  light  scattering  measure¬ 
ments,  typically  derived  by  ratioing  scattering  and  extinction,  the  optical  diameter  Dcs  is 
appropriate.  This  is 


Dea  = 


<D«> 
<  D3  > 


1/3 


\  Yli  Q<A^(mf  -  m<-i) 
ZiQt 


.jrp./ 


1/3 


(35) 


In  similar  fashion,  the  total  surface  area  of  the  particles  is  given  in  terms  of  the  Qt  by 


At 


rmir 

=  ^2  1  ?rD^(m)n<{m,t)dm 


(36) 


The  soot  volume  fraction  is  a  quantity  of  great  interest  because  it  is  usually  the  most 
important  soot  size/density  parameter  affecting  radiative  transfer.  It  corresponds  to  the 
first  moment  of  the  suspended  mass  distribution;  the  MAEROS  algorithm  gives  a  good 
representation  of  this  moment  with  a  modest  number  of  size  classes.  The  summation  over 
intrasectional  growth  coefficients  leading  to  the  time  rate  of  change  of  volume  fraction  will 
be  useful  because  it  permits  calculation  of  the  scrubbing  or  scavenging  of  gas  phase  species 
by  particle  growth.  The  analysis  described  here  is  not  suitable  for  use  at  high  pressures, 
because  continuum  effects  could  become  significant.  In  such  a  case,  growth  and  transport 
coefficimts  would  be  determined  by  interpolation,  which  would  not  seriously  impair  the 
code’s  efficiency.  The  code  could  also  readily  accommodate  a  multicomponent  particle 
analysis  (for  example,  carbon  and  hydrogen). 

2.  TVansport/Conservation  Equations  for  Particle  Sections 

The  species  conservation  equation  for  the  n^  in  the  standard  counterflow  coordinates  (Fig¬ 
ure  1)  is. 


;  I  (“«('' +''1- 


D  dni 
ne  dz 


(37) 
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where  vt  is  the  thermophoretic  velocity  and  D  the  particle  diffusion  coefficient,  respec¬ 
tively.  Equation  (371  is  multiplied  by  m  and  integrated  over  m  as  in  Equation  (19),  yielding 
(using  r  ^^(pur)  +  ^(pv)  =  0,  and  introducing  the  sectional  mass  fractions  Y<  =  pQt)' 


dYt 

^u— +,v 


3Yt 

dz 


+ 


d  /•“*  r 


I"'''’'  ■ 


dm  = 


(38) 


At  this  point,  something  must  be  known  about  the  dependences  of  vt  and  D  on  particle 
mass.  In  the  assumed  large  Knudsen  number  limit,  we  have  (Waldman,  1966;  Friedlander, 
1977;  Gomez  and  Rosner,  1994). 


-3 


.-1  /Mx  1  dT 


3  kT 

2  /2ir  R,T 

D  "V-w*- 


D' 

pVD?l^ 


(H-  T  a/8)  * 


(39) 


where  ti  is  the  gas  viscosity,  and  both  quantities  have  an  accomodation  coefficient  a  that 
is  near  unity.  It  is  seen  that  vt  is  independent  of  particle  size,  and  D  is  proportional  to 
m(~^/^).  Thus,  Equation  38  can  be  shown  to  be  equal  to  (in  similarity  form) 


"S  (l>Y,  (vt  +  VD,))  -  V  ^  +  Qf  =  0 
where  a  mass-weighted  mean  diffusion  coefficient  for  the  size  class 


(40) 


Di  = 


iD  ~  J 


and  diffusion  velocity 


VDe  = 


1 

pYr 


Dt 


B{pYi) 

dz 


(41) 


(42) 


are  empl<7ed.  The  yrtm  in  Equation  40  is  the  time  rate  of  change  due  to  aerosol 
dynamical  processes,  and  can  be  expressed  in  terms  of  the  Yr  as 
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(43) 


Q' = i  E  E  -  "'^'E 

\  i 

-j  ‘’’/3«Y?  -  „’Y<^  <">0„Yi+/.  <‘>G<Y< 

i 

+P  ^*^G4+toY<+<n  -  P  +  6ei  S(r) 


where  the  growth  and  coagulation  coefficients  are  as  derived  in  the  previous  section. 

F.  Scavenging  of  Gas  Phase  Species  by  Particulates 

Soot  growth  and  oxidation  will  deplete  local  concentrations  of  growth  species  like  acetylene 
and  benzene,  as  well  as  the  oxidating  species  oxygen  and  hydroxyl  radical.  It  is  not 
strictly  correct  to  assume  the  soot  growth  processes  represent  a  perturbation  on  the  gas 
phase  chemistry.  Scavenging  of  these  species  by  particle  processes  is  accounted  for  in  the 
following  way.  Because  p,  K  (Equation  33)  represents  the  rate  of  change  of  carbon  density 
due  to  surface  growth  and  oxidation,  the  quantities 

“  (ro-Ro,-Roh) 

So,  =  (  o - ^— )  0.  fv  (44) 

\Kg  -  ROa  -  RoH  / 

=  (rg-Ro"-Roh) 

represent  sources  and  sinks  of  solid  carbon. 

For  the  Harris- Weiner  growth  rate  assumed  in  the  sample  calculations  to  follow,  the  chem¬ 
ical  equation  can  be  represented  by 


C3H3  -J-  soot  —*■  (soot  +  20(5))  -f-  H3  (45) 

and  for  the  oxidative  processes 


O3  -I-  soot  — ►  (soot  —  C(s)]  +  O  -H  CO 


OH  -I-  soot  -♦  (soot  —  C(s))  -t-  H  -h  CO 


(46) 
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Stoichiometric  coefficients  ...  )  can  be  assigned  to  the  participating  species,  and 

molar  rates  of  change  calculated  from 

WkwJ  =  Sg  +  So,  +  SoH  (47) 

In  a  similar  way,  an  overall  chemical  equation  can  be  written  for  the  assumed  inception 
process,  and  gaseous  species  depletion  rates  calculated  from  the  inception  rate  S-,  (Equation 
15). 

III.  Computational  Results 

Sample  calculations  were  performed  for  room  temperature,  opposed  jets  of  undiluted 
methane  and  air  at  one  atmosphere.  Plug  flow  boundary  conditions  (Kee,  et  al.,  1988)  are 
employed  for  velocity  profiles  in  these  calculations.  These  boundary  conditions  provide  a 
reasonable  approximation  for  some  experimental  flame  conditions  (Ghelliah,  et.  al.,  1990). 
Table  II  gives  the  base  case  flow  parameters,  and  default  values  of  other  growth  param¬ 
eters.  Fifteen  sections  were  chosen  for  base  calculations  because  accurate  sc\>t  number 
densities  and  average  sizes  were  desired;  it  will  be  shown,  however,  that  accurate  total  vol¬ 
ume  fractions  can  be  achieved  with  far  fewer  sections.  The  aim  of  these  sample  calculations 
is  to  demonstrate  the  capabilities  of  the  program,  and  to  undertake  systematic  variations 
that  yield  more  understanding  of  the  complex  interactions  of  parameters  important  to 
soot  formation.  There  is  still  considerable  uncertainty  surrounding  the  proper  treatment 
of  inception  and  surface  growth,  and  the  focus  at  this  time  is  on  qualitative  comparisons. 

Table  II 

Base  Case  Parameters 

Vedge  =  -0.0105  gm/cm*/sec  (Yn2  =  -77,  Y02  =  -23) 

Vw»ii  =  -0208  gm/cm’/sec  (Ych4  =  1-0) 

T«dg«  —  Twaii  —  298  K 
a  =  17.2  sec"* 

Tmax  =  2000  K 

M  =  15 
PR  =  1.5 

Surface  Growth  =  HW 


Calculated  gas  phase  species  profiles  for  the  base  case  are  shown  in  Figure  2.  The  major 
growth  and  oxidative  species  are  shown,  as  well  as  species  like  CO  and  H3  that  can  be 
formed  from  growth  and  oxidation.  Oxygen  diffuses  to  the  fuel  side  of  the  fiame,  and 
is  present  at  levels  sufficient  to  affect  soot  growth,  as  will  be  discussed.  The  effect  of 
increasing  strain  rate  by  roughly  a  factor  of  five  is  shown  in  Figure  3.  While  acetylene 
is  not  strongly  affected,  the  shorter  residence  times  lead  to  reduced  benzene,  as  seen.  In 
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the  absence  of  thermal  radiation,  peak  temperatures  would  decrease  monotonicaliy  with 
increasing  strain  rate.  With  radiation  included,  peak  temperatures  are  relatively  flat  over 
the  range  of  strain  rates  considered  here  (13  to  78  see"*).  The  higher  spatial  gradients  in 
the  high  strain  rate  case  lead  to  much  higher  oxygen  concentrations. 

The  calculated  evolution  of  the  total  volume  fraction  and  the  profiles  of  the  sections  are 
shown  in  Figure  4.  The  average  size  of  the  sections  is  shown  in  nanometers.  Evolution 
toward  larger  mass  particles  occurs  due  to  surface  growth  and  coalescence.  The  low  sections 
appear  in  the  vicinity  of  the  acetylene  peak,  but  never  acquire  an  appreciable  fraction  of  the 
soot  mass.  Most  of  the  soot  growth  occurs  in  a  narrow  range  between  about  0.5  to  0.8  in  the 
normalized  flow  coordinate.  Falling  temperature,  acetylene,  and  benzene  concentrations 
give  rise  to  rapid  reductions  in  surface  growth  rates  as  the  stagnation  plane  is  approached. 
As  the  particles  grow,  convection  and  thermophoresis  push  them  in  the  direction  of  the 
stagnation  plane.  At  the  stagnation  plane,  thermophoresis  causes  some  spilling  of  large 
particles  to  the  fuel  side  of  the  plane,  but  convection  pushes  them  back.  A  balance  between 
convection  and  thermophoresis  leads  to  a  peak  on  the  fuel  side  of  the  stagnation  plane.  The 
population  of  bins  above  M  =  10  is  negligible  and  is  not  shown.  Cases  with  high  volume, 
fractions  present  in  the  largest  bin  should  be  rerun  with  an  increased  setting  for  the  upper 
limit  of  the  soot  size  range.  For  extremely  low  values  of  strain  rate  where  residence  times 
are  long  and  coalescence  has  much  time  to  promote  panicles  to  large  sizes,  adjustment 
sometimes  has  to  be  made.  As  will  be  discussed.  Figure  4  represents  exactly  the  form  of 
the  soot  volume  fraction  distribution  that  is  observed  in  experiments. 

Figure  5  shows  the  evolution  of  the  soot  mass  distribution,  with  particle  mass  as  the  in¬ 
dependent  variable,  with  distance  from  the  flame.  As  the  distribution  evolves  toward  the 
stagnation  plane,  it  takes  on  the  form  of  a  self-  preserving  distribution,  but  always  has 
so  ae  degree  of  bimodality,  with  a  secondary  peak  arising  from  inception.  A  self-preserving 
form  is  indicative  of  a  coalescence-  dominated  dynamical  process  (Friedlander,  1977).  Co¬ 
alescence  plays  an  important  role  in  determining  soot  yield  because,  although  it  is  mass- 
conserving,  it  is  destructive  of  surface  area,  and  thus  reduces  the  contributions  of  surface 
growth.  Previous  calculations  (Colket,  et.  al.,  1991)  indicated  that  these  distributions 
were  close  to  self-preserving  form  for  0.12  -  1.0  atmosphere  premixed  flames. 

Figure  6  displays  the  velocity,  velocity  gradient,  and  temperature  profiles  for  the  base 
case.  The  strain  rate,  following  Chelliah,  et.  al.  (1990),  has  been  taken  from  the  one-half 
velocity  gradient  at  the  stagnation  plane.  These  flames  are  relatively  thick,  affording  the 
long  residence  times  needed  for  soot  formation  and  significant  radiation  effects. 

Division  of  the  source  and  sink  terms  for  growth  and  oxidation  by  the  total  soot  surface  area 
leads  to  specific  rates  for  these  processes  in  units  of  cm/sec,  shown  in  Figure  7.  The  Harris- 
Weiner  (HW),  Frenklach-Wang  (FW),  and  modified  Frenklach-Wang  (MODFW)  surface 
growth  rates  have  been  discussed  previously.  When  these  are  compared  with  oxidation, 
it  is  apparent  that  only  the  Harris- Weiner  expression  results  in  significant  net  surface 
growth.  The  MODFW  expression  for  growth  is  never  in  excess  of  oxidation,  and  the 
FW  expression  (a  =  l)is  only  marginally  so  in  a  restricted  region.  Near  the  flame  front 
(peak  temperature  regime),  the  FW  mechanism  predicts  significant  growth  rates  because 
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of  increased  temperature  and  the  rise  in  concentrations  of  H-atoms.  However,  oxidation 
rates  (by  OH  and  O2)  by  far  dominate  this  growth  expression.  Calculations  including 
Soret  diffusion  did  not  alter  the  H-atom  profile  enough  to  change  any  of  these  conclusions. 

Inclusion  of  radiation  from  gas  bands  and  soot  depresses  peak  temperature  by  about  70 
K  for  the  base  strain  rate.  This  is  a  cooling  effect  of  about  3.5%,  and  probably  repre¬ 
sents  about  7%  conversion  of  flame  enthalpy  to  radiation,  the  variation  of  the  latter  with 
temperature  being  nonlinear  (Vranos  and  Hall,  1993).  This  is  large  enough  to  perturb 
significantly  the  soot  growth  species  benzene  and  acetylene,  as  shown  in  Figure  8.  The 
consequence  of  this  is  a  marked  sensitivity  of  volume  fraction  to  radiation,  as  seen  in  Fig¬ 
ure  9.  Without  radiation,  the  predicted  volume  fraction  is  nearly  a  factor  of  three  higher. 
The  higher  soot  production  with  higher  temperatures  in  a  diffusion  flame  is  consistent 
with  the  phenomenological  understanding  provided  by  Classman  (1988).  This  arises  not 
just  from  the  increased  levels  of  acetylene  and  benzene  shown  in  Figure  8,  but  also  from 
the  Arrhenius  terms  in  the  inception  and  surface  growth  expressions.  Radiation  effects 
should  decrease  monotonically  with  increasing  strain  rate,  as  the  flames  become  thinner 
and  there  is  less  residence  time  for  radiation.  This  very  strong  sensitivity  of  radiation  on 
soot  predictions  is  consistent  with  results  of  Sivathanu  and  Gore  (1993).  Thus  accurate 
calculations  of  radiation  phenomena  is  a  key  ingredient  in  the  prediction  of  soot  profiles  in 
flames.  Most  of  the  radiative  gas  cooling  is  actually  due  in  this  case  to  gas  bands,  because 
the  soot  radiation  temperature  is  relatively  low.  In  cases  with  higher  soot  volume  fraction, 
soot  radiation  would  play  a  greater  role  (see,  for  example,  Figure  16). 

In  the  radiative  power  density  formulation,  Equation  7,  the  soot  radiative  constant  C  and 
the  integrated  gas  band  intensities  (particularly  that  of  the  4.3  micron  CO2  band)  were 
made  consistent  with  Grosshandler  (1993).  The  long-wavelength  H2O  and  CO2  bands 
were  ignored.  The  low  strain  rate  flames  are  thick  enough  (Figure  6)  that  self-absorption 
effects  might  be  significant  (Hall,  1994).  The  analysis  of  such  optical  thickness  effects  is 
computationally  intensive,  however,  and  has  been  deferred  to  future  studies. 

Figure  9  also  shows  the  importance  of  oxidation.  Suppressing  it  entirely  results  in  the 
increase  in  volume  fraction  shown.  Oxidation  in  the  soot  growth  region  is  mainly  due  to 
molecular  oxygen.  Suppressing  thermophoresis  results  in  increased  soot,  as  also  shown. 
Thermophoresis  lowers  the  residence  time  in  the  growth  region,  and  thus  reduces  soot  yield. 
Without  thermophoresis,  there  is  no  spilling  of  particles  to  the  fuel  side  of  the  stagnation 
plane.  While  thermophoresis  plays  an  important  role,  it  appears  that  ordinary  diffusion 
does  not;  setting  the  sectional  diffusivities  equal  to  zero  results  in  very  little  change. 

Figures  10  and  11  show  the  analogous  effects  on  particle  number  density  and  average  size. 
Because  the  rate  of  inception  is  strongly  affected  by  temperature,  removing  radiative  loss 
results  in  much  higher  levels  of  particles,  at  least  towards  the  inception  region.  Near 
the  stagnation  plane,  there  is  less  difference  in  number  density,  presumably  due  to  a  loss 
of  sensitivity  to  initial  conditions  in  a  coalescence-dominated  regime.  Number  density 
is  decreasing  toward  the  stagnation  plane,  of  course,  because  of  coalescence.  Figure  11 
shows  interesting  effects  on  the  average  particle  size,  which  is  the  optical  diameter  Dea 
discussed  earlier.  More  rapid  surface  growth  resulting  from  removal  of  radiation,  and 
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faster  coalescence  from  higher  number  densities  results  in  much  larger  particles,  as  seen. 
The  particles  that  do  penetrate  to  the  fuel  side  of  the  stagnation  plane  evidently  do  not 
further  coalesce;  the  number  densities  and  collision  frequencies  are  too  low  for  this  to 
proceed  on  a  reasonable  time  scale.  Without  thermophoresis,  the  particles  have  more  time 
to  grow  to  a  larger  size  at  the  stagnation  plane,  as  seen. 

Variation  of  strain  rate  leads  to  the  set  of  curves  shown  in  Figure  12.  Plotted  versus 
the  reduced  coordinate,  the  total  volume  fractions  seem  to  have  the  form  of  a  similarity 
solution.  It  has  been  mentioned  that  the  peak  temperatures  are  relatively  flat  with  strain 
rate  because  of  the  effects  of  radiation.  Decreasing  benzene  and  increasing  oxygen  with 
increasing  strain  rate  are  factors  in  the  fall-off  of  volume  fraction,  but  the  main  factor  is 
probably  residence  time.  There  is  simply  less  time  for  surface  growth  to  occur  at  high 
strain  rate.  The  combined  effect  is  to  produce  an  order  of  magnitude  less  soot  with  a 
variation  of  strain  rate  of  about  six. 

Qualitative  comparisons  with  the  data  of  Zhang,  Atreya,  and  Lee  (1992)  are  quite  satis¬ 
factory,  most  notably  in  the  way  the  model  reproduces  the  profile  of  the  volume  fraction 
profile,  and  the  accumulation  of  soot  at  the  stagnation  plane.  Quantitative  comparisons 
have  not  been  attempted  because  non-zero  temperature  gradients  at  the  boundary  in  the 
experiment  is  inconsistent  with  the  modeling  assumptions. 

Because  the  surface  growth  rate  is  uncertain,  a  set  of  variations  in  which  the  nominal  HW 
rate  was  multiplied  by  factors  of  2*',  where  n  takes  on  positive  and  negative  values.  This  is 
seen  in  Figure  14.  For  very  large  enhancements  of  the  base  rate,  the  predicted  soot  yield 
begins  to  saturate  because  the  conversion  of  fuel  carbon  to  soot  becomes  large.  This  may 
be  a  trivial  effect,  refiecting  just  the  depletion  of  the  available  carbon  pool.  There  is  much 
more  relative  sensitivity  at  lower  rates.  This  is  further  exhibited  in  Figure  15.  If  scrubbing 
of  acetylene  is  removed,  the  soot  volume  fraction  grows  by  a  large  factor.  In  the  large 
surface  growth  limit,  the  predicted  volume  fraction  also  loses  sensitivity  to  inception  rate, 
as  shown  by  the  curve  in  which  the  base  inception  has  been  divided  by  the  same  factor.  In 
this  case,  the  number  density  is  markedly  reduced  in  the  inception  region,  but,  together 
with  average  size,  is  not  much  different  at  the  stagnation  plane.  The  latter  effect  possibly 
is  also  attributable  to  a  loss  of  sensitivity  to  initial  conditions  in  a  coalescence-dominated 
limit.  The  loss  of  sensitivity  to  inception  rate  in  high  growth  systems  is  consistent  with  the 
result  of  Kennedy,  et  al.  (1990).  Lower  surface  growth  rates  show  much  more  sensitivity  to 
inception  rate;  since,  in  such  cases,  predicted  soot  volume  fractions  are  nearly  proportional 
to  inception  rates. 

The  effect  on  temperature  of  many  of  these  effects  that  have  been  discussed  is  shown  in 
Figure  16.  The  elevation  of  peak  temperature  by  removal  of  radiation  discussed  earlier 
is  shown.  Enhancement  of  base  surface  growth  results  in  further  temperature  reduction 
because  soot  formation  and  thus  soot  radiation  are  enhanced.  Repeating  the  enhanced 
surface  growth  calculation  without  the  gas  phase  scrubbing  results  in  further  temperature 
reduction  because  of  even  higher  soot  levels  and  higher  radiation.  The  cooling  in  the  no¬ 
scrubbing  case  is  mitigated  to  some  degree  by  removal  of  the  exothermic  surface  growth 
process. 
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Enhanced  surface  growth  with  scrubbing  of  gas  phase  species  leads  to  the  significant  re¬ 
duction  in  the  concentrations  of  benzene  and  acetylene  shown  in  Figure  17.  Removal  of 
scrubbing  allows  these  species  to  restore  themselves,  although  these  ‘restored’  values  are 
less  than  what  they  would  have  been  if  temperature  suppression  caused  by  increased  ra¬ 
diation  was  not  included.  The  oxygen  profile  shows  similar  dramatic  effects.  Enhanced 
surface  growth  leads  to  enhanced  surface  area,  and  depletion  of  oxygen  relative  to  its  no 
scrubbing  profile.  Suppression  of  scrubbing  restores  oxygen  to  values  higher  than  its  base 
(non-enhanced  surface  growth)  profile,  because  of  lower  temperatures  due  to  enhanced  soot 
radiation.  Comparisons  in  Figures  IS,  16  and  17  provide  mple  of  the  complex 

interaction  and  interdependency  of  gas-phase  species,  particuid.Les,  bulk  temperature,  and 
radiation.  Solutions  of  one  parameter,  without  taking  into  account  effects  due  to  other 
processes,  can  be  in  error  significantly. 

The  sectional  representation  is  highly  efficient  at  predicting  soot  volume  fraction,  usually 
the  soot  parameter  of  most  interest,  with  just  a  few  sections.  Figure  18(a)  shows  how  the 
base  case  volume  fraction  profile  varies  with  number  of  sections.  For  M  >  3,  there  is  no 
significant  sensitivity  to  M!  (Strictly  speaking,  M  =  2  may  not  be  a  good  choice  since  one 
would  probably  be  violating  the  requirement  that  significant  soot  mass  not  accumulate 
in  the  highest  section).  The  solutions  for  M  >  5  are  virtually  indistinguishable.  Since 
most  of  the  soot  mass  is  due  to  surface  growth,  this  means  that  total  soot  surface  area 
is  being  accurately  represented  by  just  a  few  sections,  as  well.  The  extra  calculational 
burden  imposed  by  M  sectional  equations  is  not  significant,  and  the  sectional  coefficient 
calculations  do  not  impose  a  significant  burden  regardless  of  the  value  of  M.  Run  time 
to  convergence  has  more  to  do  with  the  accuracy  of  the  starting  solution  than  with  the 
number  of  sections.  Accurate  values  of  average  particle  size  or  number  density  cannot  be 
expected  from  so  few  sections,  however,  as  shown  for  the  average  size  in  Figure  18(b). 
Fifteen  to  twenty  sections  are  required  to  converge  on  size  or  number  density.  This  result 
is  consistent  with  studies  on  premixed  flames  (Colket  and  Hall,  1994). 


IV.  Conclusions/Sununary 

A  new  tool  for  studying  soot  formation  in  opposed  jet  diffusion  flames  has  been  developed 
by  coupling  the  dynamical  sectional  equations  for  spheroid  growth  with  the  conservation 
equations  and  complex  chemistry  of  counterflow  diffusion  flames.  The  particle  kinetics 
include  inception,  surface  growth  and  oxidation,  and  coalescence  within  the  framework  of 
the  discrete  size  spectrum,  or  sectional  representation.  The  results  of  this  analysis  indicate 
that  the  particle  growth  dynamics  and  transport,  gas  phase  chemistry,  flame  radiation 
and  flow  hydrodynamics  are  all  strongly  coupled.  Particle  inception  and  surface  growth 
oxidation  are  calculated  from  local  values  of  the  major  growth/oxidation  species  benzene, 
acetylene,  oxygen,  and  hydroxyl  radical  which  have  been  calculated  from  complex  chemical 
mechanisms.  Scrubbing  of  gaseous  species  by  the  particles  is  included  in  the  species  and 
energy  conservation  equations,  to  account  for  effects  of  soot  growth  and  oxidation.  Pcu'ticle 
transport  by  thermophoresis  and  diffusion  with  an  effective  sectional  diffusion  coefficient 
have  been  included.  Non-adiabatic  radiative  loss  by  gas  and  particle  emission  is  also 
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coupled  to  the  flow  energy  equation.  Qualitatively,  predictions  are  in  good  agreement 
with  observations  of  soot  growth  in  these  flames;  soot  is  predicted  to  appear  first  on  the 
fuel  rich  side  of  the  flame  near  the  benzene/acetyiene  peak,  and  to  accumulate  in  the 
vicinity  of  the  stagnation  plane.  Application  of  the  program  to  methane-air  flames  has 
shown  a  strong  coupling  of  thermal  radiation  and  soot  growth  as  gas  cooling  reduces  the 
concentrations  of  benzene  and  acetylene,  and  affects  Arrhenius  factors  in  the  inception 
and  surface  growth  expressions.  For  strong  surface  growth,  predicted  soot  parameters 
become  somewhat  insensitive  to  inception.  Thermophoresis  is  shown  to  reduce  local  soot 
yields  by  reducing  particle  residence  time  in  high  surface  growth  regions.  For  high  surface 
growth,  depletion  of  gas  phase  growth  species  is  shown  to  be  important  for  quantitative 
prediction  of  soot  growth.  For  nominal  surface  growth  rates,  oxidation  is  shown  to  play 
a  non-negligible  role.  The  algorithm  is  shown  to  be  highly  efficient  for  predictions  of  soot 
volume  fraction,  with  accurate  solutions  obtained  for  as  few  as  two  or  three  size  classes. 
Accurate  predictions  of  average  particle  size  or  number  density  require  larger  numbers  of 
sections.  The  work  marks  the  first  time  that  the  sectional  growth  equations  widely  used 
in  aerosol  science  have  been  coupled  to  the  solution  algorithms  for  counterflow  diffusion 
flames  that  are  important  in  combustion  science.  The  code  will  provide  a  framework  for 
evaluation  of  proposed  surface  growth  and  inception  models  as  well-posed  experimental 
data  become  available. 
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Appendix  A.  Boundary  Value  Problem  Solution  Algorithm 

Complete  specification  of  the  problem  posed  by  Equations  1-7  requires  that  boundary 
conditions  be  imposed  at  each  end  of  the  computational  domain  0<z<L.  At  z=0  (the  fuel 
jet),  we  have 


V  =  pfVf  (Al) 

U  =  0  (A2) 

Yk=Ykf,  k  =  l,2,...,K  (A3) 

Yk  =0,  k  =  K-l-l,K-l-2,...,K  +  M  (A4) 

T  =  Tf  (A5) 

and  at  z=L  (the  oxidizer  jet),  we  have 

V  =  Pox'^ox  (A6) 

U  =  0  (A7) 

Yk  =  Ykoxi  k  — 1,2, ...,K  (A8) 

Yk  =  0,  k  =  K-l-l,K-l-2,...,K  +  M  (A9) 

T  =  Tox  (AlO) 


where  Tf  and  Tox  are  the  specified  temperatures  in  the  fuel  and  ozidizer  streams;  Yk,  and 
Yk.,  are  the  specified  mass  fractions  in  the  fuel  and  oxidizer  streams  and  pfVf  =  Vedge  and 
Pox  Vox  =  Vwxii  are  the  mass  ejection  rates  from  the  fuel  and  oxidizer  jets.  By  specifying 
the  value  of  V  at  the  two  jet  surfaces  the  problem  is  overspecified  since  the  continuity 
equation  is  first-order  in  space.  To  maintain  a  well  posed  problem  the  equation 


is  introduced  and  J  solved  for  as  an  eigenvalue.  It  is  worthwhile  to  note  that  in  the  above 
formulation  the  strsdn  rate  does  not  appear  explicitly  in  the  governing  equations.  The 
magnitude  of  the  reduced  pressure  gradient  J  is  related  to  the  “stretch”  in  the  flame  due 
to  the  imposed  flow.  As  the  magnitude  of  J  increases,  so  does  the  strain  rate.  As  a 
result,  in  the  disctission  that  follows,  the  two  quantities  are  used  interchangeably.  Strictly 
speaking,  the  strain  rate  could  be  computed  directly  from  derivatives  of  the  axial  or  radial 
velocity  fields. 

Eqcstions  1-7,  Al-All  form  a  nonlinear  two-point  boundary  value  prob  "he  solution 
procedure  employs  a  combination  of  time-integration  and  an  adaptive  i  jte  difference 
method  to  obtain  profiles  for  the  dependent  variables.  The  method  has  been  discussed 
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in  detail  elsewhere  (see,  e.g.,  Smooke  (1982)  and  Smooke  et  al.  (1983,1986))  and  only 
the  essential  features  are  outlined  here.  The  goal  is  to  obtain  a  discrete  solution  of  the 
governing  equations  on  the  mesh 

X  =  (0  =  zo  <  zi  <  • .  •  <  Zm  =  L}  (-^12) 


When  the  continuous  differential  operators  are  replaced  by  finite  difference  expressions,  the 
problem  of  finding  an  analytic  solution  of  the  governing  equations  is  converted  into  one 
of  finding  an  approximation  to  this  solution  at  eeurh  point  of  the  mesh  X.  Upon  writing 
the  discrete  equations  in  residual  form,  we  seek  the  solution  U*  of  the  nonlinear  system  of 
difference  equations 

F(U)  =  0  (A13) 

Assuming  we  can  obtain  an  initial  solution  estimate  that  is  sufficiently  “close”  to  U*, 
the  system  of  equations  in  (A13)  can  be  solved  by  Newton’s  method.  We  write 

J(U‘‘)  (U'‘+^  -  U**)  = -AkF(U‘‘),  k  =  0,l,...  (A14) 

where  denotes  the  k*^  solution  iterate,  the  k^^  damping  parameter  (Deuflhard, 
1974)  (0  <  A  <  1)  and  J(U’‘)  =  3F(U‘‘)/dU  the  Jacobian  matrix.  A  system  of  linear 
block  tridiagonaJ  equations  must  be  solved  at  each  iteration  for  corrections  to  the  previous 
solution  vector.  As  we  found  in  the  solution  of  burner-stabilized  and  freely  propagating 
premixed  laminar  flames,  the  cost  of  forming  (we  use  a  numerical  Jacobian)  and  factoring 
the  Jacobian  matrix  can  be  a  significant  part  of  the  cost  of  the  total  calculation.  In  such 
problems,  a  modified  Newton  method  is  applied  in  which  the  Jacobian  is  reevaluated  pe¬ 
riodically.  An  error  estimate  that  enables  one  to  determined  when  to  form  a  new  Jacobian 
is  contained  in  Smooke  (1983). 


The  solution  of  combustion  problems,  such  as  the  counterflow  diffusion  flame,  requires  that 
the  computational  mesh  be  determined  adaptively.  The  mesh  X  is  equally  distributed 
on  the  interval  \Q,L]  with  respect  to  the  non- negative  function  W  and  the  constant  C. 
Specifically,  the  mesh  is  formed  by  employing  a  weight  function  that  equidistributes  the 
difference  in  the  components  of  the  discrete  solution  and  its  gradient  between  adjacent 
mesh  points.  Upon  denoting  the  vector  ofN  =  K-i-4  +  M  dependent  solution  components 
by  U  =  (Ui,  Ua, . . . , UnP»  we  seek  a  mesh  X  such  that 


< 


«1 


max  Ui  — 

0<i<L 


min  U 

0<»<L 


0,1,..., m  “  1 
i  =  1,2,. ..,N 


(A15) 


and 


.  ^  ,  dUi  .  dUi, 

dz  <  Tf  max  — - mm  -r— 

~  'o<i<L  dz  o<*<L  dz 


j  =  l,2,...,m-  1 
i=  1,2,...,N 


(A16) 


where  6  and  7  are  small  numbers  less  than  one  and  the  maximum  and  minimum  values  of  U{ 
and  dUi/dz  are  obtained  from  a  converged  numerical  solution  on  a  previously  determined 
mesh. 
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1  Determination  of  a  “good”  initial  estimate  for  can  be  extremely  difficult  for  counterfiow 

problems.  The  difficulty  is  due  to  the  exponential  dependence  of  the  chemistry  terms  on  the 
I  temperature  and  to  the  nonlinear  coupling  between  the  fluid  and  thermochemistry  solution 

fields.  To  alleviate  this  situation,  we  utilize  a  flame  sheet  starting  procedure  in  which 
the  equations  of  mass  and  momentum  are  coupled  with  a  Shvab-Zeldovich  equation  to 
provide  estimates  for  the  mass  flux  in  the  transverse  direction,  the  similarity  function,  the 
temperature  and  the  stable  major  gas  phase  species  in  the  flame  (Keyes  and  Smooke,  1987). 
Gaussian  profiles  are  tised  for  the  minor  gas  phase  species  and  the  soot  size  classes.  This 
starting  estimate  is  then  integrated  in  time  until  the  solution  lies  in  the  convergence  domain 
of  Newton’s  method  on  the  initial  grid.  The  size  of  the  time  steps  is  chosen  adaptively 
by  monitoring  the  local  truncation  error  of  the  time  discretization.  Time  integration  and 
Newton’s  method  are  also  employed  on  subsequently  refined  grids. 
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Coordinate  System  for  Opposed  Jet  Diffusion  Flame. 

Species  Profiles  for  Base  Case  (amuc  =  17.2  sec(-^)). 

Species  Profiles  for  Highest  Strain  Rate  Case,  (amax  =  78  sec^”^^). 

Spatial  Evolution  of  Size  Classes. 

Evolution  of  Soot  Mass  Distribution  (spatial  locations  given  in  table). 

Velocity,  Strain  Rate,  and  Temperature  Profiles  for  Base  Case. 

Comparison  of  Surface  Growth  Models,  with  Oxidation  Rate  for  Base  Case. 
Radiation  Effects  on  Growth  and  Oxidation  Species. 

Effect  of  Radiation,  Oxidation,  and  Thermophoresis  on  Soot  Volume  Fractions 
Effect  of  Radiation,  Oxidation,  and  Thermophoresis  on  Particle  Number  Density. 
Effect  of  Radiation,  Oxidation,  and  Thermophoresis  on  Average  Particle  Size. 
Variation  of  Volume  Fraction  with  Strain  Rate. 

Experimental  Soot  Distributions  and  Number  Densities  in  Methane  Opposed  jet 
Flames.  Taken  from  Figure  6  of  Zhang,  Atr^a,  and  Lee  (1992).  The  different  symbols 
show  the  effects  of  diluents.  Zn  =  z',^=fv. 

Effects  of  Surface  Growth  Augmentation.  NX  indicates  multiplication  of  base  surface 
growth  rate  by  factor  N,  inception  unchanged. 

Effects  of  Scaven^g  and  Inception  Sensitivity  at  High  Surface  Growth.  NX/M  de¬ 
notes  simultaneous  multiplication  of  base  surface  growth  by  factor  N,  and  division  of 
base  inception  rate  by  factor  M. 

Effect  of  Radiation,  Enhanced  Soot  Growth,  and  Scavenging  on  Temperature  Profile. 
Effect  of  Enhanced  Surface  Growth  and  Scavenging  on  Gas  Phase  Species. 
Sensitivity  of  Volume  Fraction  (a)  and  Average  Particle  Size  (b)  to  Number  of  Sections. 


B-33 


-0.4  -0.2 


(SSVIO  3ZIS  Nl  NOLLOVyj  3nrnoA)oo’i 


I 


B-38 


Temperature 


(m)  3dniva3dMiL 


1 - \ - 1 - “T - 1 - I 


oooooooo 

xt  ro  CN  V 

'  I  I 

(oes/ujo)  AJJ00T3A 


B-39 


B-41 


UJO)Z 


-0.2 


BdniVdldHlL 


B-49 


UTRC94-28 


Appendix  C. 

Radiative  Dissipation  in  Planar  Gas-Soot  Mixtures 
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RADIATIVE  DISSIPATION  IN  PLANAR  GAS-SOOT 

MIXTURES 
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United  Technologies  Research  Center,  East  Hartford,  CT  06108,  U.S.A. 

(Recewed  J2  May  1993) 

Abstract — A  theoretical  expression  for  the  gas  cooling  rate  (radiative  dissipation)  in  planar 
mixtures  of  soot  and  combustion  product  gases  is  derived  that  is  valid  for  all  degrees  of  optical 
thickness.  The  absorption  terms  are  reduced  to  numerical  quadrature  over  the  coordinate 
normal  to  the  flame  structure,  and  describe  soot  and  gas  self-absorption,  as  well  as  soot-gas 
spectral  overlap  effects.  Numerical  examples  illustrating  the  effect  of  radiative  cooling  on  NO 
production  in  a  model,  sooting,  opposed  jet  diffusion  flame  are  given. 


INTRODUCTION 

The  coupling  of  flame  radiation  and  flow  is  a  subject  of  renewed  interest,  with  an  acceptance  of 
the  fact  that  reacting  flow  models  are  not  complete  without  radiation  effects.  The  interest  goes 
beyond  the  ability  to  predict  heat  transfer  to  enclosure  walls.  It  is  known,  for  example,  that 
substantial  fractions  of  flame  energy  can  be  converted  to  radiation,'  and  that  the  gas  cooling 
resulting  from  this  non-adiabatic  loss  can  affect  flame  chemistry,  leading  to  inaccuracies  in 
prediction  of  pollutants  like  nitric  oxide.  Neglect  of  radiation  in  flame  hydrodynamic  simulations 
can  also  lead  to  large  errors  in  flow  velocity  and  flame  length.  These  effects  are  expected  to  be  most 
extreme  in  sooting  flames,  but  both  gas  and  soot  radiation  have  to  be  included  in  the  most  general 
model.  Quantitative  assessment  of  these  effects  requires  inclusion  of  an  energy  sink  term  in  the  flow 
energy  equation.  This  sink  term,  the  divergence  of  the  net  radiative  flux,  is  the  difference  of  an 
emission  term  valid  in  the  optically  thin  limit,  and  a  self-absorption  term  that  becomes  important 
when  optical  thickness  effects  need  to  be  accounted  for.  The  expression  for  the  emission  term  is 
generally  trivial  even  when  both  soot  and  gas  band  radiation  are  present,  and  does  not  depend  on 
flame  geometry  because  only  local  properties  are  required.  The  self-absorption  term  is  more  difficult 
because  a  solution  for  the  radiative  flux  is  implicitly  required,  and  integrations  over  all  frequencies 
and  flux  propagation  directions  are  required.  This  term  is  complicated  further  when  both  soot  and 
gas  radiators  are  present  because  of  absorption  profile  overlaps;  the  broadband  soot  absorption 
and  emission  profile  will  overlap  those  of  the  relatively  narrowband  molecular  resonances.  In  this 
paper,  a  complete  solution  for  the  radiative  dissipation  is  given  fqr  the  cne-dimensional, 
plane-parallel,  or  boundary  layer  flow  problem.  Both  gas  band  and  soot  radiation  are  considered, 
and  the  resulting  solution  is  valid  for  all  degrees  of  optical  thickness.  The  theory  will  be  illustrated 
by  example  calculations  for  a  model,  high  pressure,  sooting  counterflow  diffusion  flame.  Other 
applications  might  include  fire  radiation  from  smoke  layers,^  and  two-dimensional,  turbulent-jet 
diffusion  flames  in  the  boundary  layer  approximation.’ 

ANALYSIS 

Certain  aspects  of  this  problem  have  been  treated  by  other  investigators.^  The  main  features 
of  this  analysis  are  as  follows.  The  absorption  terms  are  derived  by  taking  hemispherical  averages 
of  frequency-dependent  solutions  of  the  equation  of  radiative  transfer,  and  integrating  them  over 
the  al^rption  bandshapes  of  both  the  molecules  and  the  soot.  The  soot  absorption  coefficient  is 
assumed  to  have  the  Rayleigh  form,  with  scattering  ignored.  For  nonuniform  paths,  the  soot  terms 
can  be  derived  without  approximation;  for  the  gas  band  terms,  a  form  of  the  mean  properties 
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assumption  is  employed.  The  gas  absorption  profiles  are  assumed  to  have  the  wideband  model 
form,  and  the  transmissivities  are  appropriate  to  the  high  pressure-broadening  limit.  More 
specifically, 

1.  It  is  assumed  that  the  molecular  absorption  bandshapes  arc  described  by  wideband  models 
for  which  the  absorption  profile  is  given  by’ 


/  X  “vP.  /-coict) -0)1;'!^ 


(1) 


where  a,j  is  the  integrated  band  intensity  of  the  jth  resonance  of  radiating  species  /  (H2O,  CO2, 
CO);  Pi  is  its  mass  density;  is  the  bandwidth  parameter  of  the  resonance,  Co  =1  for  a  band 
with  an  upper  or  lower  wavenumber  head  at  and  Co  =  2  for  a  symmetrical  band  having  a  center 
at  cul®’.  It  will  further  be  assumed  that  the  absorption  profiles  vary  much  more  rapidly  with 
frequency  than  other  characteristic  functions  such  as  Planck  function  or  soot  absorption  coefficient, 
so  that 


K^j((oY(p})  dcu  =  Co  j*  djo)  - 

=  OLi,Pif(Oif>). 


(2) 


The  molecular  transmissivity  t  is  taken  to  have  the  form  appropriate  to  the  high  pressure¬ 
broadening  limit  for  which  the  integrated  absorptance  >4;y  =  Jdto(l  -t,^)  has  the  form* 

K)  +  InK)  +  7,  (3) 

where 


Jvfii 

Ao) 


y  being  the  optical  path  length,  the  exponential  integral,  and  y  the  Euler-Mascheroni  constant. 
Its  derivative  with  respect  to  y  is  thus 


=  A©,y(l  -  e-’O/.y. 


(4) 


The  suitability  of  this  approximation  for  combustion  radiation  calculations  has  been  discussed 
previously.*  It  has  been  found  to  give  good  agreement  in  most  cases  with  more  elaborate 
calculations  based  on  narrowband  models  even  at  one  atmosphere.  Many  of  the  interesting 
applications  will  be  for  high  pressure  systems,  and  the  approximation  should  be  particularly  good 
for  these  systems.  For  nonhomogeneous  optical  paths,  evaluation  of  the  band  parameters,  which 
appear  in  Eqs.  (3)  and  (4),  are  based  on  simple,  path-averaged  temperature  and  density.*’  Overlaps 
between  molecular  bands  are  ignored,  but  the  overlaps  between  the  relatively  narrow  molecular 
resonances  and  the  broadband  soot  absorption-emission  profile  are  included.  It  will  further  be 
assumed  that  the  transmissivity  at  any  frequency  will  be  given  by  the  product  of  the  gas  and  soot 
transmissivities,  e.g.  T(aj)  =  T,(oo)r,j(co). 

2.  The  soot  absorption  coefficient  is  taken  to  have  the  form  appropriate  to  Rayleigh  spheroids 


K„{Q))  =  c,cof„  (5) 

where/,  is  the  particulate  volume  fraction,  eo  the  frequency  in  units  cm~',  and  the  constant  c,  is 
taken  to  be  independent  of  frequency.  This  is  equivalent  to  neglecting  dispersion  of  the  soot  index 
of  refraction.  The  Rayleigh  form  of  ffie  absorption  coefficient  is  suitable  to  first  order  for  aggregate 
chains  of  spheroids  as  well.'®  Light  scattering  by  soot  is  assumed  to  be  small. 

The  generalized  form  of  the  radiative  source  term  is  given  by  the  double  integral  over  wavelength 
and  direction 


<inK{<o)[I^(co)  -  m,  cu)]  =  Q^i„  - 


(6) 
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where  /fc  is  the  Planck  function,  /(fttu)  is  radiation  intensity,  and  K,(a})  =  K„(u))  +  I.,j 
Consistent  with  the  above  assumptions,  the  first  term  in  Eq.  (6)  (emission)  can  be  simply 
expressed  as"  '^ 


4ji  X  «.yP(4/  +  V,  j 

.)  Jo 


u*du 

t--\ 


cj.r, 


(7) 


where  the  factorization  assumption,  Eq.  (2),  has  been  used  for  the  molecular  resonances  ij,  the 
soot  index  dispersion  has  been  ignored  as  mentioned,  and  c,  and  Ci  are  constants  in  the  Planck 
function 

4(y, <o)  =  (exp(^)  -  >)  • 

Thus 

4i  =  c,  ((o!®>)^(exp(cja>!®Vr)  -  I )- '  =  c, (tuf )'  f  exp{ - nc.oftT).  (9) 

n  «  I 

In  the  absorptive  analysis,  a  typical  molecular  resonance  ij  will  be  assumed  to  simplify  the 
notation,  and  a  sum  over  all  ij  performed  at  the  end.  The  radiation  field  solution  is  required  as 
a  function  of  frequency  and  flux  propagation  direction  to  calculate  the  self-absorption  part. 
Referring  to  Fig.  1,  the  sum  of  the  intensities  with  direction  cosines  /i  and  -/r  with  respect  to  the 
axis  y  normal  to  the  flame  structure  can  be  formally  represented  as 


/bCK',<a)- 


y,y')<iy'  f* 

'Vp)  m  J, 


h(y\(o) 


-dT(y',y)dy' 

d(y7M)  M 
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where 


T(>',>’')  =  exp 

- 

*  K(y\oi)^^ 
y  J 

t(>’',J')  =  cxp 

- 

r^’  dv"i 

K(y\w)-^ 

y  J 

The  absorption  part  will  thus  have  terms  proportional  to  K^,,x,j{dxJd{y'lfi)), 

^t5Ti(dT,>7d(;'7;^)).  and  K„Xij{dxJd(y'lti)).  These  represent,  respectively,  gas  absorption  of  gas 
emission,  gas  absorption  of  soot  emission,  soot  absorption  of  gas  emission,  and  soot  absorption 
of  soot  emission.  Thus  +  C,  +  •  Since  the  gas  and  soot  transmissivities  have  the 

form,  respectively, 

Tj^(w)  =  exp(-^„y(ft))l>'  -y'l/|/)  =  exp(-az//i),  (11) 


and 


T,(to )  =  exp( -c.a)Jy(y;y')\y  -y'l/fx)  =  exp(-b/fi),  (12) 


where 


a  =  a,j 


Acdij 


\y-y'l 


z  =  Ztj  =  exp(-Co|tt)  -  0}f\IA0),j), 
b  =  bij  =  c,(o'^7^(y,y')(y -y% 


7Ay,y')  =  iy--yr'^[My'')^y'', 


all  of  the  terms  contributing  to  0.t»  will  be  proportional  to 

exp(-(az  +b)/ti).  (13) 

The  overbars  in  Eqs.  (11)  and  (12)  denote  path-averaged  quantities.*-’  The  contribution  from  all 
directions  of  propagation  is  given  in  terms  of  the  hemispheric  integrals 

I  dfll/(m,n)!  =  2n  J'd#i  [Eq.  10].  (14) 

Thus,  the  exponential  integral  £,  will  be  a  factor  in  all  terms  from  the  relation 

2n  I  exp( — (az  +  b)/ft)  —  =  2nE,  (az  +  b).  (15) 

Jo  fi 

The  frequency  integrations  are  generally  carried  out  using  J"  dcu-r(A(u/co)Jodz/z. 

Starting  from  the  assumed  forms  for  the  absorption  coefficients  and  transmissivities,  and  making 
use  of  the  foregoing  assumptions,  the  following  expressions  for  the  terms  contributing  to  the 
absorption  part  are  derived. 

I  ,<..»)  (16) 

=  J"^d//,(;^')W>'')/b(y(>'')£/(^  -  Ua,*)  (17) 

G„*2nc./,(y)i;<>  f*  dy''Ki:i,,{y')Iy,j{y')\y-yr'[E,{b)-E,{a-^b)]  (18) 

ij  J  -  00 
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where  a  and  b  «  b^j  are  as  given  in  Eqs.  (t  I)  and  (12),  £,  is  the  exponential  integral  of  second 
order,  and 

Qu-2itclcJ,(y)  I  (iy%(y')  f  dcuto’  £  e’'^^“'^>£,(az((u)  +  - y'l) 

J  -00  Jo  •“  I 

=  2nclcJ,{y)-^  f  dy'f,(y')T*(y')  £ 

V^2 /  J  - 00  mm  I 

-2jitrJcL/‘,()-)X(td5yy  [  iy'Hy')^,jiy')  £ 

V  J  -  QD  «'  •  I 


(E,(a+b)-E,(b)) 


P.(y')  =  nCilT(y') 

Q{y')  =  cMy-y'\- 

and  both  require  the  integral 

U(p,a,b)  =  j' dzz^Etiaz-yb).  (20) 

Series  expansions  for  this  integral  are  given  in  the  Appendix.  The  two  terms  for  Q^, 
representing  the  pure  soot  absorption  plus  a  correction  for  gas  attenuation,  respectively,  are 
derived  by  setting 

£,(«  +  6)  =  £,(b)  +  [£,(az  +  6)  -  £,(*)]  (21) 

and  assuming  that 

I  d<u(£,  (oz  +  6)  -  £, (6)]  ^  X  («  +  *)-  E,  (*)),  (22) 

J  ij 

where,  again,  a  and  6  are  as  in  Eqs.  (11)  and  (12).  Aside  from  the  factorization  assumption 
of  Eq.  (2),  Eq.  (22)  is  the  only  mathematical  approximation  made.  It  is  possible  to  derive  an 
exact  series  representation  of  this  term  for  b>a,  which  confirms  the  validity  of  this 
approximation. 

Thus  all  terms  contributing  to  field  absorption  are  expressed  in  semi>analytic  form  in  terms  of 
single  integrals  over  the  coordinate  normal  to  the  flame  structure.  These  running  integrals  are, 
however,  a  computational  burden.  The  economy  of  the  calculation  is  improved  by  restricting 
attention  to  the  strongest  molecular  lines  (HjO  6.3  and  2.7  micron,  COj  4.3  and  2.7  micron,  CO 
2.3S  micron),  which  account  for  most  of  the  radiation.  Limiting  the  terms  in  the  Planck  function 
expansion,  Eq.  (9),  to  n  =  2  is  usually  an  adequate  approximation  at  high  temperature  as  well.  This 
type  of  analysis  does  not  lend  itself  to  inclusion  of  non-black  wall  boundary  conditions.  For  black 
walls,  however,  additional  terms  in  the  expression  for  the  absorption  can  be  calculated.  As  an 
example  for  a  single  wall  at  temperature  the  following  term  is  added  to  the  absorption 
expression  for  the  pure  soot  case. 

2nc,cJAy)4\  t  -  5QJPI  ■  [^ln(l  -i-  P^QJ  -  i  ^  ’  <23) 

where 


P^  =  ncj/Tw 

Q^  =  cMy-y,\. 
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The  net  radiative  flux  to  the  wall  is,  for  pure  soot, 

^*  =  2hc.c,  j*  dy'  f  doxu^  f]  e  -y.|)/,(v') 

J-OD  Jo  *-» 

‘  Jl I'' ’  i  {(SiTrielTP  - 

K-*F;-ii(E5S¥fe)']l 

where 


(24) 


P.(y')  =  ncjny'), 

Q(y')  =  cj,\y' - yj. 

Corrections  to  Eqs.  (23)  and  (24)  for  attenuation  by  gas  band  absorption  would  be  derived  by 
proceeding  as  in  Eqs.  (21)  and  (22),  but  with  £,  replaced  by  E^. 

The  gas  band  contributions  to  net  radiative  wall  flux  may  be  expressed  as 

=  2n  I  r  dy'W)Aii}.,(y')  ly.  ->•'(-'[£,(*)  -  +  (>)1.  (25) 

Ly  J  -00 


SAMPLE  CALCULATIONS 

Sample  calculations  have  been  performed  for  a  model  opposed  jet  diffusion  flame  as  shown 
in  Fig.  1.  The  configuration  consists  of  two  opposed  fuel  and  oxidizer  jets,  leading  to  a  one¬ 
dimensional  structure  in  a  region  near  the  stagnation  streamline  where  application  of  this 
one-dimensionai  radiation  analysis  is  appropriate.  Methane  fuel  has  been  assumed,  and  the  flame 
properties  were  calculated  using  a  widely  used  opposed  jet  solver.  The  flame  was  modified  by 
adding  the  expression  for  the  radiative  source  term  derived  here  to  the  one-dimensional  energy 
equation.  The  clean  flame  calculations  are  described  in  more  detail  in  Ref.  14,  where  they  formed 
the  basis  for  a  flamelet  approach  to  predicting  NO^  concentrations  in  an  aircraft  gas  turbine 
combustor.  The  fuel  temperature  was  adjusted  to  give  a  peak  flame  temperature  approximating 
that  of  jet  fuel,  and  the  pressure  of  lO.Satm  is  representative  of  such  a  combustor  at  the  cruise 
condition.  In  order  to  assess  the  possible  role  of  soot  in  the  radiative  transfer,  a  sooting  region 


Fig.  2.  Adiabatic  temperature  profile  in  model  opposed  jet  and  assumed  location  of  sooting  region.  Strain 

rate*  20  sec*'. 
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has  been  synthetically  added  to  the  clean  flame  solution  as  shown  in  Figs.  I  and  2.  Soot  will 
tend  to  form  on  the  fuel  side  of  the  flame,  and  wtll  build  up  at  the  stagnation  plane  due  to 
thermophoretic  transport.'*  “  Calculation  cf  the  actual  soot  profile,  of  course,  requires  a  soot 
kinetics  and  transport  model,  but  this  was  net  available  for  this  study.  Thus,  a  simple  picture  of 
the  soot  profile  is  assumed  in  which  it  has  been  assumed  to  be  uniformly  distributed  between  the 
stagnation  plane  and  a  point  on  the  fuel  side  where  nucleation  might  be  expected  to  be  intense. 
Following  Ref.  17,  we  assume  that  nucleation  will  occur  at  an  elemental  carbon  mixture  fraction 
of  0.1,  where  the  adiabatic  temperature  is  about  I67S  K.  This  results  in  the  soot  profile  located 
as  shown  in  Fig.  2;  the  effective  radiation  temperature  of  the  soot  is  about  1550  K,  nearly  900  K 
lower  than  the  peak  gas  temperature.  The  average  soot  volume  fraction  is  then  treated  para¬ 
metrically.  The  synthetic  nature  of  this  calculation  must  be  emphasized;  with  methane  fuel  at  this 
temperature  the  flame,  in  reality,  would  be  expected  to  be  relatively  clean.  The  intent  here  is  to 
gain  some  idea  of  the  relative  importance  between  soot  and  gas  band  radiation  for  reasonable 
ranges  of  the  soot  volume  fraction.  The  soot  is  coupled  to  the  gas  flame  solution  only  through 
temperature  reduction  resulting  from  the  radiation;  in  reality,  the  gas  cooling  due  to  radiation 
could  affect  the  pyrolysis,  nucleation,  and  surface  growth  processes  that  grow  the  soot.  An  accurate 
treatment  of  this  problem  requires  that  temperature-dependent  soot  growth  processes  be  added  to 
the  chemical  kinetics  mec  tism.  The  base  gas  kinetics  mechanism  used  here,  where  the  effect  of 
radiation  on  NO  concentration  is  of  interest,  is  a  variation  of  that  proposed  by  Glarborg  et  al." 
The  213  reactions  and  57  species  of  that  mechanism  have  been  culled  to  126  and  30,  respectively. 
The  data  set  is  limited  to  Q  chemistry,  and  NOj  kinetics  have  been  ignored. 

Opposed  jet  solutions  are  characterized  by  the  strain  rate  (denoted  here  by  the  symbol  a),  which 
is  a  measure  of  the  velocity  gradient  normal  to  the  flame  structure.  It  is  an  indication  of  how  fast 
reaction  products  are  pulled  away  in  the  direction  parallel  to  the  flame  sheet.  High  values  of  strain 
rate  correspond  to  short  gas  residence  times,  thinner  flame  structures,  and  somewhat  lower  peak 
temperatures,  leading  to  reduced  importance  of  radiation.  Low  values  of  strain  rate  are  character¬ 
ized  by  the  opposite  attributes,  and  it  will  consequently  be  these  that  are  examined  here.  In 
particular,  a  base  case  value  of  20  sec'*  is  used  for  this  study;  the  adiabatic  temperature  profile 
is  that  shown  in  Fig.  2.  For  low  values  of  strain  rate,  chemical  equilibrium  radiating  species  profiles 
will  be  a  reasonable  approximation,  and  lower  values  of  strain  rate  can  be  simulated  by  stretching 
the  coordinates  in  proportion  to  a  ~  Figures  3  and  4  are  based  on  fixing  the  adiabatic  solution 
for  a  =  20sec~',  stretching  the  coordinates  in  this  way,  and  calculating  the  radiative  dissipation 
according  to  Eqs.  (7)  and  (16>-<I9).  Figure  3  shows  how  calculated  absorption  and  emission  vary 
with  strain  rate  for  a  clean  flame  condition.  It  can  be  seen  that  optical  thickness  effects  are 
important  even  for  the  a  =  20  sec  case,  where  the  characteristic  flame  thickness  is  about  1-2  mm. 
Figure  4  shows  the  analogous  results  for  a  sooting  flame  where  a  soot  volume  fraction  of  I0*~** 
has  been  assumed.  At  this  value  of  volume  fraction,  the  soot  contributions  start  to  dominate  the 
radiation  process  at  this  pressure. 


Fig.  .i.  Variation  of  radiative  dissipation  profiles  with  strain  rate  for  non-soottng  case. 
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Displacement  (cm)  Displacement  (cm) 

Fig.  4.  Variation  of  radiative  dissipation  profiles  with  strain  rate  for  soot  loading  of  10*'’’. 


Inclusion  of  Eqs.  (7)  and  (16)-(19)  as  a  power  sink  term  in  the  energy  equation  of  the  opposed 
jet  solver  leads  to  the  predicted  variation  of  radiative  loss  fraction  with  soot  volume  fraction  shown 
in  Fig.  S.  As  stated,  soot  radiation  starts  to  dominate  for  volume  fractions  in  excess  of  about 
with  predicted  loss  fractions  rising  to  about  20%  for  the  highest  values  of  volume  fraction  that 
could  be  reasonably  expected.  Soot  self-absorption  becomes  important  at  the  highest  volume 
fraction.  Radiative  loss  fraction  is  here  defined  as  the  volume  integrated  net  radiative  power  density 
divided  by  the  volume  integrated  chemical  enthalpy  release  rate.  These  curves  are,  of  course,  highly 
sensitive  to  the  assumed  thickness  and  effective  radiation  temperature  of  the  soot  zone.  Calculations 
at  other  values  of  strain  rate  are  qualitatively  similar,  with  the  overall  radiative  loss  fraction 
decreasing  with  increasing  strain  rate. 

The  gas  cooling  effects  resulting  from  the  radiative  loss  are  shown  in  Fig.  6.  The  percentage 
decrease  in  flame  temperature  will  generally  be  much  less  than  the  radiative  loss  percentage  because 
the  heat  capacity  is  nonlinear,  and  because  reductions  in  temperature  will  be  offset  to  some  extent 
by  recombination  of  radicals.  In  the  case  of  soot  radiation,  the  soot  will  be  radiating  from  a  lower 
temperature  region,  and  the  effect  on  the  temperature  there  is  shown  in  Fig.  6  for  the  maximum 
value  of  volume  fraction  considered.  Cooling  of  the  flame  zone  will  then  occur  by  thermal 


Fig.  S.  Variation  of  predicted  radiative  loss  fraction  with  soot  loading  for  optically  thin  and  thickness-cor¬ 
rected  cases. 
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Fig.  6.  Temperature  profiles  for  adiabatic,  gas  band  radiation,  and  gas  band  plus  soot  radiation. 


Displacement  (cm) 

Fig.  7.  Predicted  effect  of  radiative  cooling  on  NO  profiles. 


conduction;  the  percentage  reduction  in  peak  flame  temperature  is  far  less  than  the  20%  radiative 
loss  fraction  for  this  case.  As  stated,  however,  the  temperature  reductions  in  the  soot  zone  are 
probably  large  enough  to  affect  the  soot  growth  processes,  leading  to  a  complicated  coupling 
between  radiation  and  growth. 

The  effect  of  these  temperature  reductions  on  the  profiles  of  NO  is  shown  in  Fig.  7.  The  profile 
resulting  from  the  soot  radiative  losses  can  be  contrasted  with  that  resulting  from  assuming  a  local 
non-adiabatic  loss  that  is  everywhere  20%  of  the  local  enthalpy  release  rate.  Because  the  main  effect 
of  the  soot  is  to  cool  regions  where  the  NO  is  not  at  its  peak  levels,  soot  radiation  has  a  much 
smaller  effect  on  NO  reduction  in  these  flames. 
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APPENDIX 


The  integral 

(•I  (*1  (*00  •(-“  + 

£/(p,  a,6)  =  j  dzz^Ei(az  +  b)=  \  dt — j— 

has  the  following  series  representations: 


Azz^Eyiftz  -I-  h)  =  £,(a  -F ft) /(p  -I- 1) 


i: 

rc-(.>siy/  ^  Yg,-i(g+») 

1  i\a+b}  (p-f-l),+  , 

e-'r(p -H)y*(p -f  l,a)  ft  a  V  ««(<=>+ ft) 

(P  +  l)  '(P  +  l)\ft +«/  I  +  ^/  (P  +  U,.,! 


(p>-l) 

(ft  >a) 
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where 


(Al) 


(A2) 


(“)«+!  =  f(M-)-n  -I-  !)/£(«) 

«,(«)  =  n!(l  -Fir  -I-  u^/2]  -F  ■  •  u"ln\) 


and  r  and  y*  are  the  gamma  function  and  a  form  of  the  usual  incomplete  gamma  function,  y, 
respectively. 

The  first  series  is  derived  by  performing  the  z -integration  in  Eq.  (Al)  first,  using” 


i: 


dzz‘’e  =  i4,(p, af,  1)  = 


y(p  -F  l.g  -F  1) 


(A3) 


y  is  now  expanded  in  a  power  series,®  and  the  r -integral  is  done  term-by-term.  The  second  series 
is  derived  by  using  Eq.  (A3)  and  integration  by  parts,  followed  by  term-by-term  integration 
over  t. 
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Appendix  D. 

Modeling  of  Soot  Aerosol  Formation  in  Opposed  Jet  Diffusion  Flame 

(abstract) 


To  be  submitted  for  presentation  at  Fourth  International  Aerosol  Conference,  Los 
Angeles,  August  28-September  2, 1994. 


MODELING  OF  SOOT  AEROSOL  FORMATION  IN  OPPOSED  JET 
DIFFUSION  FLAMES.  R.  J.  Hall.  United  Technologies  Research  Center,  East 
Hartford,  CT  06108,  M.  D.  Smoote,  Department  of  Mechanical  Engineering, 
Yale  University,  New  Haven,  CT  06250,  and  M.  B.  Colket,  United  Technologies 
Research  Center,  East  Hartford,  CT  06108. 

Sectional  dynamical  equations  for  spheroid  growth  have  been  coupled  to 
conservation  equations  for  opposed  jet  diffusion  flames  with  complex  chemistry 
and  transport  to  provide  an  a^ytic id  tool  fOT  studies  of  soot  formation.  The 
ao’osol  dynamical  ^nations  use  the  sectional  method  of  Gelbard  and  Seinfeld 
(1980),  with  inception,  surface  growth/oxidation,  and  coalescence.  Particle 
transport  by  thermophoresis  and  diffusion  is  included.  The  free  molecule  regime  is 
aq)propriate  to  atmospheric  pressure  flames,  affording  simplification  of  the 
dynamical  and  trans^rt  coefficients.  The  particulate  species  conservation 
equations  are  coupl^  to  the  hydrodynamic  and  gas-phase  species  conservation 
equations  of  the  opposed  jet  or  counterflow  solver  of  Smool^,  Puri,  and  Se^dri 
(1986).  Particle  inception  is  based  on  an  expression  involving  calculated  local 
concentrations  of  acetylene  and  benzene  derived  from  complex  chemical 
mechanisms.  Surface  growth  occurs  at  a  rate  that  is  governed  by  temperature  and 
local  concentration  of  ^  major  growth  species  acetylene,  and  likewise  oxidation  is 
governed  by  concentrations  of  oxygen  and  hydroxyl  radical  Thennal  radiation  by 
particles  and  gas  phase  species  is  represented  by  a  sink  term  in  the  hydrodynamic 
energy  equation.  Scrubbing  of  gas  jAase  species  by  particle  growth  and  p^cle 
thermochernistry  are  also  treated,  llie  residt  is  a  strongly  coupled  gas-particulate 
system.  This  work  marks  the  first  time  that  the  sectional  growth  equations  widely 
u^  in  aerosol  science  have  been  coupled  to  the  solution  schemes  for  counterflow 
diffusion  flames  that  are  important  in  combustion  science. 

Comparisons  of  model  predictions  with  experimental  data  are  presented.  Both 
theory  and  experiment  show  soot  inception  to  occur  on  the  fuel  side  of  the  flame, 
with  a  monotonic  increase  of  volume  fraction  or  suspended  aerosol  mass  as 
thermophoresis  and  convection  push  the  particles  back  in  the  direction  of  the  fuel 
stream.  Maximum  volume  fraction  occurs  theoretically  at  the  stagnation  plane,  with 
few  particles  able  to  cross  to  the  fuel  side,  as  observed  experimentally.  Extensive 
parametric  variations  exhibiting  the  importance  of  thermal  radiation,  particulate 
scrubbing  of  gaseous  species,  thermophoresis,  and  oxidation  will  be  presented. 

Gelbard,  F.  and  Seinfeld,  J.  H.  (1980).  J.  Coll.  Int  Sci.  78:485-501. 

Smooke,  M.  D.,  Puri,  I.  K.  and  Seshadri,  K.  (1986).  Twenty-First  Symposium 
(International)  on  Combustion,  The  Combustion  Institute,  p.  1783-1792. 
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Appendix  E. 

Predictions  of  Soot  Foraiation  in  Counterflow  Diffusion  Flames 

(abstract) 


To  be  submitted  for  presentation  at  Twenty-Fifth  International  Symposium  on 
Combustion,  Irvine,  CA,  July  31-August  5, 1994. 


PREDICTIONS  OF  SOOT  FORMATION  IN  COUNTERFLOW  DIFFUSION 
FLAMES.  M.  D.  Smoc^e,  Department  of  Mechanical  Engineering,  Yale  University, 
New  Haven,  CT  06250,  R.  J.  Hall,  Uni^  Technologies  Research  Center,  East 
Hartford,  CT  06108,  and  M.  B.  Colket,  United  Technologies  Research  Center,  East 
Hartford,  CT  06108. 

Dynamical  equations  for  spheroid  growth  have  been  coupled  to  conservation  equations 
for  opposed  jet  diffusion  flames  with  complex  chemistry  and  transport  to  provide  an 
anal^^cal  tool  for  studies  of  soot  formation.  The  aerosol  dynamic^  equations  use  the 
discrete  size  spectrum  or  sectional  method  of  Gelbard  and  Seinfeld  (1980),  with 
inception,  surface  growth/oxidation,  and  coalescence.  Particle  transport  by 
thermophoresis  and  diffusion  is  included  in  the  sectional  transport  equation.  The 
analysis  is  carried  out  in  the  free  molecule  regime  appropriate  to  atmospheric  pressure 
flames,  affording  simplification  of  the  dynamical  and  transport  coefficients.  The 
particulate  species  conservation  equations  are  coupled  to  the  hydrodynamic  and  gas- 
phase  species  conservation  equations  in  elliptic  form  (Smooke  and  Giovangigli,  1992). 
Particle  inception  is  based  on  an  expression  involving  calculated  local  concentrations  of 
acetylene  and  benzene  derived  from  complex  chemical  mechanisms.  Various  surface 
grow^  models  have  been  evaluated,  and  oxidation  by  oxygen  and  hydroxyl  radical  is 
included.  Optically  thin  thermal  radiation  by  particles  and  gas  phase  species  is 
represented  by  a  sink  term  in  the  hydrodynamic  energy  equation.  Scrubbing  of  gas 
phase  species  by  particle  growth  and  particle  thermochemistry  are  also  treat^.  The 
result  is  a  strongly  coupl^  gas-particulate  system  in  which  there  can  be  significant 
perturbation  of  the  gas  phase  chemistry  by  particulate  growth  and  radiation. 

The  model  predictions,  which  are  in  good  qualitative  agreement  with  methane  flame 
experiments,  show  the  importance  of  thermal  radiation,  particulate  scrubbing  of 
gaseous  species,  and  thermophoresis.  Parametric  variations  of  surface  gro\^  and 
inception  rates  wiU  be  presented  which  provide  information  on  the  sensitivity  of  soot 
yield  to  these  growth  parameters.  The  sectional  algorithm  is  shown  to  be  highly 
efficient,  yielding  accurate  soot  volume  fractions  with  as  few  as  three  sections. 

Gelbard,  F.  and  Seinfeld,  J.  H.  (1980).  J.  Coll.  Int.  Sci.  78:485-501. 

Smooke,  M.  D.,  and  Giovangigli,  V.  (1992).  Twenty-Fourth  Symposium 
(International)  on  Combustion,  The  Combustion  Institute,  p.  161. 
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A  detailed  chemical  kinetic  mechanism  for  the  pyrolysis  of  toluene  and  the 
formation  of  polyaromatic  hydrocarbons  has  'been  developed.  This  chemical  kinetic 
model  is  consistent  with  shock  tube  data  with  mass  spectral  identification  of 
intermediate  species.  This  kinetic  analysis  provides  a  better  understanding  of 
how  aromatic  rings  decompose  and  how  polyaromatic  species  grow.  In  addition,  a 
previously  developed  soot  formation  code,  based  on  a  sectional  aerosol  model  for 
predicting  soot  inception,  growth,  and  oxidation  in  a  premixed  flame  has  been 
fully  integrated  into  an  opposed-jet,  diffusion  flame  code.  The  new  code 
includes  effects  due  to  radiation  from  both  gaseous  species  and  particulatea  as 
well  as  scavenging  of  species  by  soot.  The  code  treats  particle  transport 
including  thermophores is.  Predictions  from  a  low  strain  rate,  lightly-sooting, 
methane- fueled,  opposed-jet,  diffusion  flame  are  included  in  this  report.  Per¬ 
turbation  studies  demonstrate  the  importance  of  fully  integrating  soot  production, 
radiation,  and  scavenging  in  order  to  reasonably  predict  bulk  parameters  (cont'd) 


soot  formation  modeling,  chemical  kinetics  of 
toluene  pyrolysis,  formation  mechanisms  and  thermodynamics  of 
polyaromatic  hydrocarbons,  soot  formation  in  opposed  jet  flames 
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Abstract  (Continued) 


such  as  temperature,  as  well  as  species  concentrations,  sooting  levels  and 
radiation  loads. 


